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FOREWORD 


This  report  was  prepared  by  Inyong  Ham,  Professor  of  Industrial 
Engineering,  Department  of  Industrial  Engineering,  Pennsylvania  State 
University,  University  Park,  Pennsylvania,  in  compliance  with  Contract 
DAAF01-70-C-1069  under  the  direction  of  the  Research  Directorate, 
Weapons  Laboratory,  with  R.  A.  Kirschbaum  as  Project  Engineer. 

This  work  was  authorized  as  part  of  the  Manufacturing  Methods  anc* 
Technology  Program  of  the  U.S.  Army  Materiel  Command,  which  is  admini¬ 
stered  by  the  U.S.  Army  Production  Equipment  Agency. 
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ABSTRACT 

The  work  described  in  this  report  was  initiated  by  personnel 
of  the  Research  Directorate,  Weapons  Laboratory,  U.S.  Army  Weapons 
Command,  for  application  of  computer  optimization  techniques  to  the 
study  of  machining  parameters. 

Two  computer  methods  for  industrial  optimization  of  machining 
condit^ns  are  described  and  demonstrated.  The  PERFORMANCE  INDEX 

METHOD  (PIM)  requires  only  shop  data  for ““Chi~”8 p—™ TION  *r  ° 
pieces  produced,  and  number  of  tool  changes.  The  PRODUCTION 
OPTIMIZATION  METHOD  (POM)  requires  tool  life,  time,  and  cost  data. 
Both  are  designed  to  refine  the  initial  data  input  with  shop  test 
data  obtained  during  normal  production,  as  related  to  one  or  mo*«  °f 
three  production  objectives:  minimum  unit  cost,  maximum  production 
rate  and  maximum  profit  rate.  The  computer  programs  are  constructed 
for  use  by  shop  personnel  with  little  knowledge  of  mathematics  or 

computers. 

Both  methods  are  rapid  and  economical,  and  the  programs  can  be 
processed  by  either  in-plant  or  remote  computer  facilities.  The  user 
is  given  all  information  needed  to  install  the  programs  and  adapt 
them  to  his  purposes. 
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PART  1.  INTRODUCTION 


Though  computer  assisted  systems  such  as  numerical  control  and  adap¬ 
tive  control  can  greatly  increase  the  efficiency  of  material  removal,  we 
still  need  to  improve  the  economics  of  this  industrial  operation.  It  is 
the  purpose  of  this  report  to  show  how  optimization  methods  employing  the 
computer  can  reduce  costs  by  finding  the  most  effective  machining  condi¬ 
tions  for  a  given  production  objective. 

NEW  METHODS  OF  OPTIMIZATION  ANALYSIS 

Described  here  are  two  practical  computer  methods  for  analyzing  op¬ 
timum  cutting  conditions,  developed  and  tested  in  the  Machinability 
Laboratory,  Department  of  Industrial  Engineering,  at  The  Pennsylvania 
State  University.  Both  are  designed  to  determine  from  actual  production 
data  the  best  combination  of  machining  parameters  in  relation  to  minimum 
cost  per  piece,  maximum  production  rate,  or  maximum  profit  rate. 

The  PERFORMANCE  INDEX  METHOD  requires  only  information  on  total  ma¬ 
chining  time,  number  of  pieces  produced,  and  number  of  tool  changes, 
which  machine  operators  can  easily  supply.  It  is  used  when  tool  life 
equations,  time  study  data,  and  cost  data  are  not  readily  available. 

This  method  is  based  on  a  concept  of  response  measured  directly  from  pro¬ 
duction  operations  (a  performance  index) ,  by  which  the  optimization  cri¬ 
teria  are  defined.  The  optimum  set  of  conditions  is  then  found  by  a 
mathematical  procedure  employing  production  test  results  from  the  shop. 

The  PRODUCTION  OPTIMIZATION  METHOD  requires  tool  life,  time,  and 
cost  data  generated  by  production  tests.  The  tool  life  data  are  analyzed 
in  relation  to  other  parameters,  with  repetitive  feedback  from  the  shop, 
until  adequate  tool  life  data  are  obtained.  Optimum  cutting  conditions 
are  determined  from  the  results  of  this  analysis. 

Input  to  the  computer  programs  for  these  methods  is  expressed  in  or¬ 
dinary  machining  terms,  and  the  functional  engineering  decisions  are  in¬ 
corporated  in  the  program  design.  After  the  basic  programs  have  been 
constructed,  they  can  therefore  be  used  with  varying  job  conditions  by 
anyone  who  is  capable  of  supplying  the  necessary  shop  data.  Both  methods 
are  economical  of  computer  time. 

ADVANTAGES  OF  COMPUTER  ANALYSIS 

The  primary  advantage  of  these  methods  is  that  they  can  easily  pro¬ 
vide  information  that  has  heretofore  been  difficult  or  impossible  to  ob¬ 
tain  in  the  manufacturing  plant.  It  is  a  further  advantage  that  the 
output  of  the  computer  programs  is  printed  instructions  and  tables,  imme¬ 
diately  ready  for  use  by  plant  personnel. 

A  production  engineer  who  attempts  optimization  analysis  by  conven¬ 
tional  means  has  complicated  work  to  do.  Though  he  may  have  abundant 
machinability  data  to  draw  upon,  from  sources  such  as  the  U.S.  Army  (14) 
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and  the  U.S.  Air  Force  (13),  those  data  are  from  experimental  investiga¬ 
tions  and  production  case  studies  performed  elsewhere.  They  are  not 
likely  to  apply  exactly  to  a  specific  shop  job  with  its  own  complex  vari¬ 
ables,  and  extending  them  to  a  particular  production  case  involves  many 
difficulties  and  precautions.  Analyzing  tool  life  data  by  conventional 
means  requires  lengthy  calculation  of  theoretical  and  empirical  relations. 
Numerous  efforts  to  simplify  the  practical  application  of  machinability 
data  with  worksheets,  nomographs,  and  tables  have  not  achieved  much  sim¬ 
plification.  At  best  the  calculated  solutions  are  approximations,  which 
the  engineer  has  to  interpret  from  his  own  experience  and  judgment.  More¬ 
over,  since  optimum  conditions  vary  from  job  to  job,  machine  to  machine, 
and  tool  to  tool,  he  has  to  figure  his  way  through  the  whole  routine  re¬ 
peatedly. 

So  it  is  not  surprising  that  in-plant  optimization  of  cutting  condi¬ 
tions  is  not  a  common  practice.  Consequently,  manufacturers  are  often 
unaware  that  inefficient  machining  may  be  a  factor  obstructing  their 
efforts  to  reduce  the  cost  of  their  products  on  a  highly  competitive  mar¬ 
ket.  And  this  inefficiency  can  persist  even  where  processes  are  computer 
controlled. 

Earlier  work  by  this  investigator  and  others  (3,6,7)  suggested  the 
value  of  the  computer  in  analyzing  optimum  cutting  conditions.  But  the 
programs  then  developed  were  limited  to  particular  applications  or  were 
based  on  generalized  machinability  ratings  or  given  tool  life  data.  They 
were  of  little  help  in  production,  and  were  not  comprehensive  enough  to 
include  some  of  the  interrelated  parameters. 

The  two  methods  presented  here  make  broad  use  of  the  computer's  abil¬ 
ity  to  integrate  many  parameters  simultaneously.  These  programs  can  cal¬ 
culate  rapidly  the  best  combination  of  machining  conditions  for  any  actual 
cutting  job,  in  terms  of  the  tooling  available  and  the  desired  production 
objective.  They  are  designed  for  effective  use  by  production  personnel 
who  know  the  practical  problems  and  by  the  machinists  who  perform  the 
cutting  operations.  In  plants  where  a  digital  computer  is  not  yet  stand¬ 
ard  equipment,  the  programs  can  be  serviced  by  remote  and  time-sharing 
systems . 

In  addition  to  its  cost-saving  potential,  computer  optimization  can 
be  invaluable  in  analyzing  proposed  machine  work.  New  materials,  new 
tools,  and  new  machinery  present  many  manufacturing  problems  that  can  be 
solved  by  this  kind  of  analysis.  It  can  take  much  of  the  guesswork  out 
of  estimating  the  physical  and  economic  feasibility  of  undertaking  a  job. 

COST  OF  COMPUTER  OPTIMIZATION 

The  following  cost  data  for  the  demonstrations  of  the  new  optimiza¬ 
tion  methods  described  herein  are  a  conservative  indication  of  the  cost 
of  obtaining  a  complete  computer  analysis  for  an  actual  production  case. 

Running  the  computer  program  for  the  PERFORMANCE  INDEX  METHOD  re¬ 
quired  3.32  seconds  of  compiler  time,  2.17  seconds  of  execution  time,  and 
a  total  time  of  9.0  seconds  including  input  and  output  times.  The  input 
was  847  cards  (program  and  data),  and  the  output  was  1285  printed  lines. 
With  a  time  cost  of  $0.25  per  second  for  9.0  seconds  and  a  printout  rate 
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of  $0,001  per  line  of  output,  the  total  cost  of  the  computer  analysis  was 
$3.55.  This  represents  the  maximum  cost  of  computation.  In  production, 
less  than  one-third  of  this  output  would  be  generated  at  each  iteration, 
with  a  corresponding  reduction  of  computer  time. 

When  the  computer  program  for  the  PRODUCTION  OPTIMIZATION  METHOD  was 
forced  to  generate  five  iterations  during  one  computer  run  in  the  simula¬ 
tion  mode,  the  total  computer  time  was  62  seconds  and  1486  lines  of  out¬ 
put  were  printed.  At  the  given  time  and  printout  rates,  the  cost  of  the 
complete  analysis  was  about  $17. 

Using  the  information  provided  in  this  report,  a  competent  FORTRAN 
programmer  could  adapt  either  of  these  programs  to  a  particular  user's 
needs  in  roughly  a  month's  time.  Such  tailoring  would  be  likely  to  yield 
significant  savings  of  computer  time,  input  data,  and  out  records,  since 
the  programs  are  now  written  to  cover  a  fairly  wide  range  of  conditions. 

The  shop  tests  called  for  in  both  methods  can  be  carried  out  as  nor¬ 
mal  production  machining.  Though  the  cutting  conditions  to  be  tested  will 
obviously  not  be  equally  efficient,  they  will  all  be  within  a  range  ap¬ 
proaching  the  optimum. 

Savings  of  production  costs  as  the  result  of  computer  optimization 
will  be  proportional  to  the  difference  between  previous  shop  practice  for 
the  analyzed  operation  and  the  computed  optimum  cutting  conditions. 

VERSATILITY  OF  THE  METHODS 

The  PERFORMANCE  INDEX  METHOD  can  easily  be  generalized  for  applica¬ 
tion  to  operations  other  than  turning,  like  milling  and  drilling,  and  also 
for  multitool  operations.  The  program  described  in  this  report  is  spe¬ 
cifically  designed  to  optimize  two  parameters,  cutting  speed  and  feed,  but 
it  can  be  expanded  to  add  depth  of  cut  or  other  related  factors. 

Since  the  PRODUCTION  OPTIMIZATION  METHOD  depends  on  analysis  of  tool 
life  data  for  a  single  cutting  tool,  it  could  not  be  extended  to  a  multi¬ 
tool  operation  without  complete  revision  of  the  mathematical  procedure. 
Wherever  sufficient  tool  life,  cost,  and  time  data  can  be  obtained,  this 
method  is  recommended  for  plant  use  because  it  requires  fewer  tests  than 
the  performance  index  method  and  tends  to  give  more  exact  results. 
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PART  2.  PERFORMANCE  INDEX  METHOD 


The  performance  index  method  (PIM)  computes  optimum  machining  condi¬ 
tions  for  a  given  production  index,  from-  test  data  limited  to  machining 
time,  the  number  of  pieces  produced  per  time  unit,  and  the- number  of  tool 
changes  during  that  time. 

PIM  is  designed  to  achieve  three  practical  objectives:  a)  to  obtain 
accurate  results  with  an  economical  amount  of  shop  testing,  b)  to  provide 
a  flexible  system  that  can  accommodate  varying  job  conditions,  and  c)  to 
make  optimization  analysis  as  "automatic"  as  possible,  so  that  it  can  be 
carried  out  by  machine  shop  personnel. 

The  PIM  computer  program  automatically  performs  the  mathematical 
operations  indicated  in  Section  2.2  and  makes  the  necessary  computations. 
Use  of  the  program  needs  no  understanding  of  the  mathematical  processes 
and  little  knowledge  of  computers.  All  of  the  input  and  output  is  in  the 
form  of  common  machining  data. 

The  structure  of  the  program  and  the  simple  language  code  for  com¬ 
municating  with  the  computer  are  discussed  in  Section  2.3.  Sections  2.4 
and  2.5  demonstrate  PIM  procedure  for  simulation  tests  and  for  actual 
production  analysis. 


2.1  CONCEPT  AND  TECHNIQUES 

The  performance  index  method  employs  the  performance  index  concept, 
as  its  name  implies,  and  the  general  methods  of  response  surface  tech¬ 
niques.  We  shall  begin  by  defining  these  basic  means.  Applying  them  to 
efficient  solution  of  the  problems  of  optimization  analysis  required  spe¬ 
cific  adaptive  techniques,  described  in  Section  2.2. 

PERFORMANCE  INDEX  CONCEPT 

PIM  derives  its  flexibility  from  the  performance  index  concept  (12) . 
A  performance  index  can  be  any  measurable  response  of  the  machining  oper¬ 
ation  to  changes  in  the  variables  or  parameters  involved  (11) . 

Since  the  purpose  of  optimization  analysis  is  to  improve  the  eco¬ 
nomics  of  machining,  the  logical  choice  of  indexes  would  be  unit  cost, 
production  rate,  or  profit  rate,  alone  or  in  some  combination.  These  are 
the  indexes  selected  for  the  PIM  program,  but  any  other  measurable  re¬ 
sponse  can  be  substituted. 

RESPONSE  SURFACE  TECHNIQUES 

A  response  surface  is  the  surface  described  by  one  or  more  variables 
and  a  measured  response  to  those  variables.  In  the  PIM  program,  the  re¬ 
sponse  surface  is  described  by  the  cutting  speed,  the  feed  rate,  and  the 
chosen  performance  index.  Other  variables  can  be  selected,  such  as  depth 
of  cut. 
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The  term  ’’response  surface  techniques"  encompasses  a  large  group  of 
methods  for  exploring  response  surfaces  to  find  some  optimum  point  (4,12). 
In  general  these  techniques  stem  from  the  assumption  that  the  response 
surface  cannot  be  described  by  an  exact  equation,  or  that  existing  data 
are  not  complete  enough  to  permit  accurate  analysis  by  statistical  methods 
such  as  regression  techniques.  The  usual  method  is  to  search  for  the  ap¬ 
proximate  optimum  point  by  running  tests  at  several  points  on  the  surface. 

The  aim  of  every  response  surface  technique  is  to  locate  the  optimum 
point  in  the  most  efficient  manner.  But  what  is  efficient  for  one  case 
may  not  be  efficient  for  another.  A  technique  employing  only  a  few  test 
points  may  find  an  approximation  that  is  acceptable  in  some  applications; 
in  a  different  problem,  an  inaccurate  guess  might  reduce  the  advantage  of 
using  so  few  test  points.  So  it  is  the  wide  variety  of  applications,  each 
with  its  own  characteristics,  that  explains  why  so  many  techniques  have 
been  devised. 

Unfortunately,  none  of  the  existing  techniques  appear  to  be  simple 
and  flexible  enough  for  use  in  production.  Porter  and  Summers  (12)  report 
methods  designed  for  an  adaptive  control  application.  They  involve  seek¬ 
ing  the  optimum  point  by  performing  a  series  of  single-point  tests,  ana¬ 
lyzing  the  results,  and  selecting  new  test  points  on  the  basis  of  those 
results.  This  technique  seems  quite  impractical  because  of  the  number  of 
tests  needed  and  the  inaccuracy  of  the  final  result.  The  "Box-Wilson" 
method,  described  by  Duncan  (4) ,  takes  fewer  tests  and  produces  sufficient 
accuracy,  but  it  apparently  requires  the  services  of  someone  capable  of 
making  statistical  decisions  as  the  testing  progresses.  The  necessity  for 
keeping  a  skilled  statistician  on  hand  greatly  restricts  the  utility  of 
this  method. 

PIM  also  uses  response  surface  techniques  to  locate  the  optimum 
point.  The  computer  program  finds  both  the  area  containing  the  optimum 
point  and  the  precise  position  of  the  optimum  point.  How  efficiently  it 
does  this  work,  regardless  of  the  placement  of  the  initial  test  points  on 
the  response  surface,  is  demonstrated  by  the  simulation  examples  in  Sec¬ 
tion  2.4.  Moreover,  PIM  maintains  its  efficiency  under  varying  job  con¬ 
ditions,  and  it  is  capable  of  optimizing  other  manufacturing  operations 
in  addition  to  metal  cutting. 


2.2  MATHEMATICAL  PROCEDURE 

Mathematical  relations  for  the  performance  index  method  are  presented 
in  six  steps,  corresponding  to  the  analytic  steps  of  the  PIM  computer  pro¬ 
gram.  The  basic  computational  technique  is  the  statistical  method  of 
least  squares  regression.  Since  the  functional  procedures  are  written 
into  the  program,  where  they  operate  automatically,  the  information  in 
this  section  will  be  of  interest  only  to  installing  engineers  and  inves¬ 
tigators  of  optimization  analysis. 

As  here  described,  the  PIM  program  has  written  into  it  a  subroutine 
for  simulation  tests.  This  subroutine,  which  was  used  for  the  simulation 
examples  in  Section  2.4,  employs  certain  mathematical  relations  that  are 
not  called  for  in  the  main  program.  They  are  discussed  at  the  end  of  this 
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section.  The  simulation  subroutine  can  also  be  useful  in  preliminary 
computer  studies  of  real  production  cases. 

STEP  1.  SELECTION  OF  INITIAL  TEST  POINTS 

The  first  step  in  analyzing  a  given  response  surface  is  to  choose 
the  initial  test  points.  Predetermined  tool  life  data  are  usually  needed 
to  make  this  choice,  but  PIM  is  intended  particularly  for  situations  where 
reliable  tool  life  data  are  not  available.  The  computer  program  itself 
therefore  recommends  the  initial  test  points.  Provision  is  made,  however, 
for  selection  of  the  test  points  by  an  experienced  engineer  when  adequate 
data  are  at  hand. 

PIM  has  two  basic  functions:  a)  to  find  the  general  area  of  the  op¬ 
timum  point,  and  b)  to  pinpoint  the  actual  position  of  that  point.  Since 
much  testing  can  be  saved  by  locating  the  initial  tests  in  the  area  where 
the  optimum  point  is  likely  to  occur,  the  program  is  designed  to  look  for 
expected  properties  of  the  response  surface. 

Certain  assumptions  are  made  concerning  the  shape  of  the  response 
surface  being  tested.  The  main  assumption  is  that  the  contours  of  equal 
values  of  the  performance  index  selected  (cost  per  piece,  production  rate, 
or  profit  rate)  will  be  parabolic.  This  means  that  the  problem  must  be 
constrained,  and  that  the  optimum  point  will  lie  on  the  constraint  that 
closes  the  open  end  of  the  parabola.  PIM  therefore  chooses  initial  test 
points  in  the  area  of  the  constraint  and  the  open  end  of  the  parabolic 
contours,  as  indicated  in  Fig.  1. 

The  computer  picks  nine  points,  geometrically  arranged.  For  a  re¬ 
gression  surface  with  two  variables,  the  number  of  points  should  not  be 
less  than  six.  But  the  regression  will  give  a  much  more  accurate  fit  if 
about  nine  points  are  supplied,  and  the  regression  equation  will  fit  the 
performance  index  surface  most  accurately  if  the  test  points  are  in  a 
geometric  pattern. 

Figure  1  shows  the  general  scheme  for  seeking  the  optimum  point. 

The  curves  are  the  expected  contours  for  constant  values  of  the  chosen 
performance  index,  cost  per  piece.  The  constraints  are  the  maximum  and 
minimum  speed  and  feed  and  the  maximum  horsepower  available  for  the  ma¬ 
chining  operation  (2,8).  The  computer  chooses  the  nine  initial  test 
points  so  that  they  fall  within  the  horsepower  constraint  and  span  the 
upper  half  of  the  usable  feed  values  and  the  middle  two-fourths  of  the 
usable  speed  values. 

As  mentioned  earlier,  the  program  user  can  select  the  nine  initial 
points  himself  if  he  has  enough  reliable  data  or  if  he  can  make  a  very 
good  educated  guess.  The  points  do  not  have  to  fall  in  the  area  chosen 
by  the  computer  program.  No  matter  where  they  are,  the  PIM  program  will 
probably  find  the  optimum  point,  but  convergence  to  the  true  optimum  may 
be  slower. 

STEP  2.  COMPUTATION  OF  THE  PERFORMANCE  INDEX 

We  have  defined  the  performance  index  as  any  measurable  response  to 
the  variables  in  an  operation.  After  the  tests  at  the  recommended  points 
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PERFORMANCE  INDEX  CONSTRAINT 
(equal  cost  per  piece  curves) 


FIG.  1.  Scheme  for  seeking  optimum  point  with 
nine  initial  test  points.  Usable  cutting  spaed 
and  feed  ranges  are  restricted  by  performance 
index  and  horsepower  constraints. 


have  been  run,  the  following  information  from  each  test  must  be  fed  back 
to  the  computer: 

Number  of  parts  produced  N 

Number  of  tool  changes  N 

Time  period  of  test,  min  T 

In  the  basic  PIM  program  the  performance  index  is  a  linear  ccmfcina- 
tion  of  the  production  rate  and  the  reciprocal  of  the  cost  per  piece. 

But  it  can  be  any  other  measurable  response  derived  from  the  test  results. 
The  only  restriction  is  that  there  will  be  just  one  optimum  point  in  .the 
region  of  testing. 

Theoretically  s,  the  performance  index  PI  can  be  defined  as 

PI  =  QPr  +  (1  -  Q)(l/Cu)  (1) 

where  Q  =  performance  index  criterion 
~  production  rate  ~  N^/T 

=  cost  per  piece 
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The  cost  per  piece  or  unit  cost  is  determined  by 

C  =  [ (RLO)T  +  (TLC)N  ]/N 
u  tc  p 

in  which  KLO  =  labor  and  overhead  rate,  $/min 
TLC  =  tooling  cost,  $/edge 

The  performance  index  will  consist  of  Q  x  100%  of  the  production  rate 
and  (1  -  Q)  x  100%  of  the  reciprocal  of  the  cost  per  piece.  The  value  of 
Q  will  usually  be  either  1  or  0,  but  it  can  be  some  intermediate  value  if 
a  specific  percentage  ratio  between  the  production  rate  and  the  cost  per 
piece  is  known. 

STEP  3.  REGRESSION  ANALYSIS 

When  all  of  the  production  data  for  the  initial  test  points  have  been 
collected s  we  are  ready  to  begin  analyzing  them  by  the  process  of  least 
squares  regression. 

The  following  set  of  points  is  stored  by  the  program: 


Vi 

F 1 

PIi 

V2 

0 

• 

Fa 

PI2 

V9 

Fg 

Pl9 

where  V  =  cutting  speed,  fpm,  or  spindle  speed,  rpm 
F  =  feed,  ipr 

To  fit  a  response  surface  to  these  points  by  multiple  regression 
techniques,  we  need  a  generalized  equation  that  will  fit  the  experimental 
points.  If  we  assume  that  the  surface  is  well  behaved  and  that  there  is 
only  one  optimum  point,  the  following  polynomial  equation  is  a  good  rep¬ 
resentation  of  the  response  surface  within  the  test  area: 

PI  =  b  +  b  V  +  b  F  +  b  V2  +  b  F2  +  b  VF  (2) 

12  3  4  5  6 

The  error  difference  e  between  the  experimental  value  y  of  PI  and 
the  theoretical  value  of  PI  given  by  the  right-hand  side  of  Eq.  1  can 
then  be  expressed  as 

e  =  y  -  (b  +  b  V  +  b  F  +  b  V2  +  b  F2  +  b  VF)  (3) 

1  2  3  4  5  6 

The  purpose  of  regression  analysis  is  to  minimize  the  sum  of  the 
squared  error.  Since  there  are  actually  n  such  equations,  one  for  each 
test  point,  Eq.  3  can  be  written  in  matrix  form: 

E  =  Y  -  XB  (A) 
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in  which 


*  > 

ei 

*  > 
yi 

b,  1 

1  Vi  Fi  Vi  Ff  ViFi 

E  = 

e2 

Y  = 

y2 

B  = 

cvi 

&  • 

X  = 

1  v2  f2  V2  Fi  v2f2 

•  ••••• 

e 

y 

• 

• 

b6 

•  ••••• 

1  V  F  V2  F2  V  F 

l  n  J 

n  J 

-  4 

n  n  n  n  n  n 

v. 

In  this  form,  we  can  express  Ee2  ,  the  criterion  to  be  minimized,  as  4>. 

Then 

<J>  =  (Y  -  XB)  '  (Y  -  XB) 

=  Y'Y  -  ZB'X'Y  +  B'X'XB 

and  to  find  the  set  of  b's  that  will  minimize  <j),  the  derivative  of  (j)  is 
taken  with  respect  to  each  b^: 

Scjj/Sb^^  =  -2X'Y  +  2X'XB  =  0 
X'XB  =  X’Y 

B  =  (X'X^X'Y 

Given  the  surface  described  by  Eq.  2,  it  is  now  possible  to  find  the 
optimum  point  on  that  surface. 

STEP  A.  OPTIMIZATION  OF  PERFORMANCE  INDEX 

The  optimum  point  could  be  located  by  a  mathematical  optimizing 
method.  But  since  we  are  dealing  with  discrete  points  and  a  high  speed 
computer,  the  most  efficient  way  to  locate  the  desired  point  is  a  simple 
exhaustive  search.  The  optimum  point  estimated  by  the  regression  surface 
is  easily  found  by  calculating,  with  Eq.  2,  the  value  of  PI  at  every  pos¬ 
sible  combination  of  feed  and  speed  within  the  known  constraints. 

Because  the  theoretical  regression  surface  as  defined  by  Eq.  2  has 
only  one  optimum  point,  that  is,  only  one  maximum  or  minimum,  the  search 
for  the  optimum  can  extend  beyond  the  region  of  the  tested  points.  If 
the  optimum  point  lies  within  the  region  of  testing,  the  regression  equa¬ 
tion  will  estimate  its  position  fairly  accurately.  If  it  lies  outside  the 
tested  area,  the  regression  equation  will  give  a  general  estimate  of  the 
position  or  will  at  least  show  the  direction  to  look  for  further  test 
points . 

STEP  5.  FURTHER  TESTS  FOR  OPTIMIZATION 

After  an  initial  estimate  of  the  optimum  point  has  been  obtained, 
further  testing  is  necessary  either  to  confirm  that  point  or  to  find  a 
new  estimate.  The  position  of  the  estimated  point  will  determine  what 
values  the  program  will  recommend  for  the  next  set  of  tests.  The  loca¬ 
tion  of  these  test  points  on  the  response  surface  will  depend  on  whether 
or  not  the  estimated  optimum  point  is  on  a  constraint. 
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Optimum  Point  on  a  Constraint.  As  stated  before,  we  expect  that  the 
optimum  point  will  be  found  on  a  constraint.  The  program  is  therefore 
structured  to  look  for  this  probability  and  take  advantage  of  it.  . 


When  the  estimated  optimum  point  occurs  on  one  of 
straints  (maximum  feed  or  maximum  horsepower) ,  further 
along  that  constraint.  A  curvilinear  regression  model 
test  points  to  the  constraint,  with  reference  again  to 
tion  for  the  surface  (Eq.  2)  : 


the  upper  con¬ 
testing  is  done 
is  used  to  fit  the 
the  general  equa- 


PI  =  b  +BV+bF+bV2+bF2+bVF 

1  2  3  4  5  6 


Where  F  =  F  , 
max 


PI  *  b  +  b  V  +  b  F  +  b  V2  +  b  F2  +bVF 

12  3  max  4  5  max  6  max 


(5) 


We  can  group  the  constants  and  redefime  them  as 


b  ’  =  b  +  b  F  +  b  F2 
i  i  3  max  5  max 


b’  =  b  +  b  F 
226  max 


b'  =  b 

3  4 


The  result  is  a  new  regression  equation: 

PI  =  b’  +  b'V  +  b'V2  (6) 

1  2  3 


If  there  is  a  horsepower  constraint,  the  maximum  horsepower  HP  is 
defined  by 


HP  =  [  (ttD)FV(DC)S,  ]  /  (ef  f ) 
m  np 


(7) 


where  D 
V 
DC 


hp 

eff 


diameter  of  workpiece,  in. 
spindle  speed,  rpm 
depth  of  cut,  in. 
specific  horsepower,  hp/in.3/min 

efficiency  of  machine  tool,  % 


We  can  rewrite  Eq.  7,  grouping  the  constants  in  a  new  constant  C: 

C  =  (HP  )  (eff)/[irD(DC)S,  ]  =  FV  ®r  C  =  FV  (8) 

m  np 

\ 

Solving  for  feed, 


F  =  C/V 


and  substituting  the  F  into  Eq.  2  results  in 

PI  =  b  .+  b  V  +  b  (C/V)  +  b  V2  +  b  (C/V)2  +  b  C  (9) 

1  2  3  4  5  6 

which  can  be  rearranged  as 

PI  =  b"  +  b"V  +  b" (1/V)  +  b"V2  +  b" (1/V2)  (10) 

1  2  3  4  5 

If  the  values  of  V  are  fairly  large,  say  in  the  hundreds  or  more,  the 
terms  with  V  or  V2  in  the  denominator  will  be  significantly  small  and  can 
be  dropped  from  the  equation,  which  then  becomes 

PI  -  b'  "  +  b"  'V  +  b'  1  'V2  (11) 

1  2  3 

Equations  6  and  11,  since  they  are  of  the  same  form,  can  be  joined 
as  one  continuous  regression  model  for  the  purpose  of  this  program.  Note, 
however,  that  Eq.  11  is  the  projection  of  the  intersection  of  the  horse¬ 
power  plane  and  the  PI  surface,  on  the  plane  formed  by  V  and  PI  coordi¬ 
nates.  It  is  not  the  actual  curve,  but  the  projection  of  the  curve  on  a 
plane.  As  such,  it  is  usable  in  the  program  because  the  projection  clear¬ 
ly  illustrates  the  maximum  value  of  PI  that  is  being  sought. 

The  test  points  must  be  taken  along  the  constraints  in  question,  but 
the  tests  will  be  fitted  to  the  regression  model  given  by  Eq.  6.  The  con¬ 
straints  that  will  be  tested  for  the  optimum  value  of  PI  are  shown  in 
Fig.  2.  The  circled  black  dot  indicates  the  estimated  position  of  the 
optimum  point,  and  the  open  circles  indicate  the  approximate  position  of 
the  points  recommended  by  the  program  for  further  testing. 

Optimum  Point  Not  on  a  Constraint.  In  this  case,  one  of  two  things 
may  have  occurred.  Either  the  final  position  of  the  optimum  has  not  yet 
been  pinpointed,  or  the  performance  index  surface  is  not  as  predicted. 

If  either  of  these  conditions  exists,  the  program  will  recommend  test 
points  closely  grouped  around  the  previous  estimate,  and  the  surface  re¬ 
gression  model  (Eq.  2)  will  again  be  used  to  fit  these  new  tests  to  the 
surface. 

STEP  6.  STOPPING  CRITERION  AND  ERROR  ANALYSIS 

The  stopping  criterion  for  the  PIM  program  is  that  whenever  the 
optimizing  search  produces  the  same  point  twice  in  succession,  the  program 
stops.  It  is  based  on  the  opinion  that  if  two  sets  of  test  points  yield 
the  same  result,  that  result  can  be  accepted  with  reasonable  safety.  Be¬ 
cause  the  program  is  dealing  with  discrete  points,  this  restriction  is 
not  so  tight  as  it  may  seem.  When  the  computer’s  exhaustive  search  lo¬ 
cates  an  optimum  point,  the  actual  optimum  probably  does  lie  somewhere 
within  the  area  close  to  the  computer's  recommendation. 

How  tight  the  stopping  criterion  really  is  depends  on  the  distance 
between  the  related  speed  and  feed  settings.  Figure  3  illustrates  this 
situation.  The  open  circles  indicate  a  few  of  the  possible  feed  and  speed 
settings,  and  the  circled  black  dot  indicates  the  recommended  optimum 
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point.  The  shaded  block  represents  the  area  that  can  be  assumed  to  con¬ 
tain  the  true  optimum. 

Steps  1  through  5  are  repeated  until  the  stopping  criterion  is  met. 
The  resulting  optimum  conditions  can  be  taken  to  represent  the  best  esti¬ 
mate  of  the  true  optimum,  within  the  limits  of  experimental  error. 

Definite  problems  will  arise  if  the  experimental  or  production  test 
data  contain  a  large  amount  of  error  or  variation.  During  production 
tests,  every  effort  should  be  made  to  reduce  the  effects  of  operating 
variables  such  as  the  operator’s  speed  and  time  out  for  repairs.  Large 
errors  will  distort  the  shape  of  the  performance  index  surface,  produce 
variations  in  the  optimum  point  from  one  run  to  the  next,  and  disrupt  the 
stability  of  the  analytic  method. 

The  amount  of  error  may  be  so  great  that  the  program  will  not  find 
an  optimum  point  according  to  the  stopping  criterion,  but  will  oscillate 
around  the  point.  If  the  area  of  oscillation  is  small,  the  user  can  de¬ 
cide  to  pick  the  optimum  at  any  point  within  it.  With  so  much  error  in 
the  input  data,  an  optimum  point  located  by  the  computer  would  not  be 
very  accurate  anyway.  To  obtain  the  accuracy  of  which  PIM  is  capable, 
errors  in  the  test  data  should  be  minimal. 

Statistical  error  might  also  affect  the  stopping  criterion.  This 
kind  of  error  relates  to  how  well  the  regression  surface  fits  the  actual 
performance  index  surface.  Least  squares  regression  can  minimize  statis¬ 
tical  error,  but  cannot  completely  eliminate  it.  If  the  area  of  testing 
is  small  enough  so  that  the  enclosed  surface  is  smooth  and  well  behaved 
and  has  only  one  optimum  point,  the  fit  of  the  regression  surface  should 
not  introduce  any  significant  error. 

COMPUTER  SIMULATION  TESTS 

The  validity  and  flexibility  of  PIM  were  proved  by  the  computer  sim¬ 
ulation  tests  described  in  Section  2.4.  These  tests  required  a  subroutine 
for  simulation  of  the  machining  data  (MACSIM) ,  to  evaluate  the  performance 
index.  Note  that  this  subroutine  is  used  only  in  simulation  tests.  It 
is  not  a  part  of  the  PIM  main  program  for  production  application,  which 
is  independent  of  predetermined  tool  life,  cost,  and  time  study  data. 

To  compute  the  performance  index,  as  previously  defined  for  an  actual 
run,  we  start  with  a  known  tool  life  equation  and  a  set  of  suitable  cost 
data  and  time  study  data.  From  this  information,  the  production  rate  and 
unit  cost  can  easily  be  calculated  by  the  following  mathematical  proce¬ 
dure. 


Let  us  take  as  the  tool  life  equation 
VTaF^  =  C 

where  V  =  cutting  speed,  fpm  a,£  =  exponents 
T  =  tool  life,  min  C  =  constant 

F  =  feed,  ipr 
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and  where 


V  -  ttDN/12  or  N  -  12V/ttD 

in  which  D  «  diameter  of  workpiece,  in. 

N  =  spindle  speed,  rpm 

We  first  solve  for  tool  life  T: 

T  =  (C/VF^)1^ 

We  next  calculate  T^,  machining  time  per  piece: 


T  =  L/NF 
c 


where  L  =  length  of  workpiece,  in. 

We  can  now  calculate  the  total  cost  per  piece  C^: 

C  =  C  T  +CT  +(CTT)/T+CT/T 
u  om  oc  octc  ec 

where  Cq  =  labor  rate  plus  overhead,  $/min 
=  machine  handling  time,  min 
T  =  tool  changing  time,  min 
=  tool  cost,  $/edge 

And  we  can  also  calculate  P  ,  the  production  rate: 


P  =  1/ [T  +  T  +  T  (T  /T)] 
r  m  c  ct  c 


The  performance  index  PI  can  then  be  computed  as  defined  in  Eq.  1: 

PI  =  QP  +  (1  -  Q) (1/C  ) 
r  u 

To  lend  realism  to  the  computer  simulation  tests,  a  "random  percent 
age  error"  was  added  to  each  value  of  PI  returned  to  the  subroutine,  by 
writing  into  the  simulation  program  an  error  noise  of  5%  or  10%. 


2.3  COMPUTER  PROGRAM  FOR  PIM 

The  computer  program  for  the  performance  index  method  implements  the 
mathematical  procedure  explained  in  Section  2.2.  It  employs  the  language 
of  FORTRAN  IV  and  can  be  translated  with  any  FORTRAN  compiler.  The  com¬ 
puter  used  in  the  developmental  work  was  the  IBM  360,  Model  67,  in  the 
Computation  Center  at  The  Pennsylvania  State  University.  Since  the  pro¬ 
gram  does  not  need  a  large  amount  of  storage,  it  can  be  modified  to  fit 
many  other  computer  systems.  It  can  also  be  adapted  to  optimize  other 
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machining  operations,  such  as  milling  and  drilling,  simply  by  changing 
the  input  and  output  formats. 

PIM  is  designed  to  provide  maximum  flexibility,  at  the  same  time  en¬ 
suring  efficient  and  accurate  performance.  As  in  any  other  analytic  proc¬ 
ess  involving  variables  and  testing,  there  is  a  possibility  that  the 
program  might  not  find  an  optimum  point  in  a  reasonable  time.  But  because 
of  the  nature  of  the  optimization  problem  and  the  analytic  techniques,  the 
probability  of  failure  should  be  minimal  within  the  practical  limits  of 
testing. 

This  section  describes  decisions  written  into  the  program,  several 
options  available  to  the  user,  the  coding  of  input  data,  and  finally  the 
structure  of  the  PIM  computer  program* 

DECISIONS  AND  OPTIONS 

The  PIM  optimization  procedure  requires  certain  engineering  deci¬ 
sions.  To  make  the  program  usable  by  shop  personnel  who  may  be  unfamiliar 
with  the  mathematical  method  and  may  know  little  about  computers,  these 
decisions  are  written  into  it  and  will  operate  automatically  unless  they 
are  changed.  For  flexibility,  several  options  can  be  exercised  to  suit 
specific  applications. 

Distance  Between  Adjacent  Test  Points.  Though  the  analysis  can  ac¬ 
commodate  either  a  stepped  or  a  stepless  machine  tool,  a  step-size  factor 
is  required  for  both  cases.  All  step  sizes  are  expressed  in  terms  of  a 
possible  number  of  speed  and  feed  settings,  not  as  differences  between 
the  actual  values  of  the  variables .  That  is,  the  distance  between  any 
two  points  is  some  number  of  settings,  usually  a  fraction  of  the  total 
number  available.  This  fraction  was  set  at  1/4  of  the  total  number,  as 
the  result  of  a  series  of  tests  that  proved  it  to  ensure  the  greatest  de¬ 
gree  of  generality  in  the  use  of  the  program.  By  this  engineering  deci¬ 
sion,  the  program  automatically  computes  step  sizes  as  1/4  of  the  possible 
number  of  step  settings,  either  actual  or  assumed. 

Constraints .  Speed,  feed,  and  horsepower  constraints  are  written 
into  the  program  to  ensure  that  the  recommended  optimum  condition  will 
always  be  within  the  capability  of  the  machine  tool  to  be  used.  Though 
an  optimization  analysis  can  be  performed  without  any  constraints,  the 
results  will  be  “unrealistic.  When  values  for  specific  horsepower  and  ma¬ 
chine  efficiency  are  unobtainable,  the  PIM  program  can  be  run  without  the 
horsepower  constraint.  On  the  other  hand,  it  is  possible  to  read  in  an 
additional  constraint  such  as  surface  finish,  if  necessary,  to  further 
restrict  the  usable  feeds  and  cutting  speeds  (9) . 

Initial  Test  Points.  PIM  is  programmed  to  recommend  the  initial  test 
points.  But  if  the  necessary  data  are  available,  the  user  can  bypass  this 
part  of  the  program  and  pick  the  initial  test  points  himself  (Section  2.5, 
Selection  of  Initial  Test  Points).  He  must  bear  in  mind  however,  that  at 
least  six  points  should  be  chosen,  and  that  they  should  be  arranged  in  a 
closely  grouped  geometric  pattern  to  secure  the  best  regression  fit  (Sec¬ 
tion  2.2,  Step  1) , 
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Substitution  of  Test  Points.  The  user  also  has  the  option  of  replac¬ 
ing  any  test  point  recommended  by  the  computer  with  a  point  of  his  own 
choosing.  When  the  recommended  points  are  all  on  the  upper  constraints, 
the  new  point  must  also  be  on  an  upper  constraint.  Any  replacement  must 
be  made  with  care  not  to  destroy  the  mathematical  procedure  and  the  gen¬ 
eral  trend  of  the  method.  Instead  of  a  substitution,  any  one  or  two 
points  may  simply  be  dropped,  as  long  as  the  remaining  points  are  not 
fewer  than  the  minimum  number  needed  for  the  regression  (six.  points  for 
surface  regression,  three  points  for  curvilinear  regression). 

Choice  of  Spindle  Speed  or  Cutting  Speed.  As  the  program  is  now 
written,  the  two  parameters  affecting  the  performance  index  are  spindle 
speed  N,  rpm,  and  feed  F,  ipr.  Spindle  speed  was  chosen,  rather  than  the 
cutting  speed  V,  fpm,  commonly  used  in  analyses  of  machining  operations, 
because  shop  personnel  are  often  more  inclined  to  think  in  terms  of  it. 
Spindle  speed  can  easily  be  converted  to  the  corresponding  cutting  speed, 
and  either  is  acceptable  in  the  program.  But  if  cutting  speed  is  chosen, 
the  horsepower  equation  in  Step  5  of  the  mathematical  procedure  (Section 
2.2)  will  have  to  be  rewritten  accordingly,  and  the  computer  must  be  in¬ 
structed  to  print  out  cutting  speed,  fpm,  instead  of  spindle  speed,  rpm, 
as  it  now  does.  The  user  has  also  the  option  of  substituting  some  other 
variable,  like  depth  of  cut,  for  either  of  the  two  parameters. 

Number  of  Parameters.  The  mathematical  procedure  can  handle  any  two 
discrete  parameters.  In  fact,  it  is  adaptable  to  any  number  of  parame¬ 
ters,  and  to  any  performance  index  that  has  only  one  maximum  value  in  the 
test  area.  For  more  than  two  parameters  the  program  would  have  to  be 
suitably  adjusted,  of  course.  But  even  three  parameters  will  create  prac¬ 
tical  problems,  because  the  number  of  test  points  needed  for  regression 
increases  greatly  with  each  additional  parameter. 

INPUT  DATA  CODE 

The  following  tabulation  lists  and  defines  the  code  used  for  input 


data  in  the 
notation  is 

PIM  program, 
also  shown. 

For  convenience,  the  corresponding  mathematical 

Input  Code 

Notation 

Definition 

IDENT 

Type  of  machine  tool: 

=  1  if  speed  is  stepless.  Computer 
lates  artificial  steps. 

=  0  if  speed  is  stepped.  Steps  are 
on  data  cards. 

calcu- 

read  in 

VLOW 

Vi 

Lowest  usable  spindle  speed,  rpm. 

VFLOW 

Vi 

Lowest  usable  cutting  speed,  fpm,  corre¬ 
sponding  to  Vi. 

VHIGH 

vn 

Highest  usable  spindle  speed,  rpm. 

VFHIGH 

V 

n 

Highest  usable  cutting  speed,  fpm, 
sponding  to  V^. 

corre- 
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Input  Code  Notation 


Definition 


FLOW 

Fi 

Lowest  usable  feed,  ipr. 

FHIGH 

F 

n 

Highest  usable  feed,  ipr. 

[Value  of  FHIGH  must  be  the  value  of  a  feed  setting.  If  the 
desired  maximum  feed  is  about  0.0100  ipr,  the  closest  feed 
setting  must  be  chosen  (say,  0.0102  ipr).  Otherwise  program 
will  not  execute  properly.  If  FHIGH  is  left  blank  or  read  in 
as  zero,  program  will  set  as  FHIGH  the  highest  feed  read  in. 
Similarly,  if  zero  values  are  read  in  for  FLOW,  VLOW,  or  VHIGH, 
computer  will  set  lowest  or  highest  value  read  in.] 


SHP 

Shp 

Specific  horsepower,  hp/in.3/min. 

HPM 

HP 

m 

Maximum  horsepower. 

EFF 

eff 

Efficiency  of  machine  tool,  percentage  in 
decimal  form  (for  example,  0.5). 

DC 

DC 

Maximum  depth  of  cut  expected,  in. 

RLO 

RLO 

Labor  and  overhead  rate,  $/min. 

TLC 

TLC 

Tool  cost,  $/edge„ 

Q 

Q 

Performance  index  criterion.  Usually  0  or 
1,  but  may  be  some  value  between. 

D 

D 

Diameter  of  workpiece,  in. 

VI 

VA 

VCH 


[The  following  three  variables  apply  only  if  IDENT  =  1.] 


V. 

i 


V 


ch 


Lowest  available  spindle  speed,  rpm. 

Highest  available  spindle  speed,  rpm. 

Desired  size  of  steps  between  V.  and  V  . 

i  a 


L 


[The  following  two  variables  apply  only  if  IDENT  =  0.] 

£  Number  of  spindle  speed  settings. 


V(I)  V± 


Values  of  spindle  settings,  rpm,  i  =  1  to  L 
Must  be  in  ascending  order,  Vi  lowest. 


M  M  Number  of  feed  settings. 

F(I)  F.  Values  of  feed  settings,  ipr,  i  =  1  to  M. 

Must  be  in  ascending  order,  Fi  lowest. 


KT 


Sequence  number  for  successive  computer 
runs.  KT  =  0  for  initial  run. 
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(CONTINUED) 


Input  Code 
IR 


Notation 


Definition 


NPT  or  IF 


VMAX  or  NOP  V 

max 


FMAX  or  FOP  F 

max 


Code  for  type  of  regression  and  use  of  test 
results.  Value  for  IR  Is  found  on  print  of 
previous  run.  IR  =  0  for  initial  run. 

Number  of  test  results  to  be  read  in.  NPT 
and  IF  are  used  interchangeably.  Value  is 
found  on  print  of  previous  run.  If  not  all 
test  points  recommended  were  run,  value 
should  be  actual  number  of  tests.  IF  =  0 
for  initial  run. 

Previously  calculated  value  of  optimum 
spindle  speed,  found  on  print  of  previous 
run.  NOP  =  0  for  initial  run. 

Previously  calculated  value  of  optimum  feed, 
found  on  print  of  previous  run.  FOP  =  0  for 
initial  run. 


[The  following  five  variables  are  read  in  only  if  NPT  ^  0.] 


VL(I) 


V. 

i 


Value  of  spindle  speed  used  in  a  test  run. 


FD(I) 


Value  of  feed  used  in  a  test  run. 


P(I) 


Number  of  parts  produced  in  time  T  at  test 
values  of  V^  and  F^. 


T  (I) 
TP  (I) 


T. 

l 


tci 


COMPUTER  PROGRAM 


Time  period  of  test  run,  min. 

Number  of  tool  changes  in  time  T  at  test 
values  of  V^  and  F^. 


The  PIM  computer  program  is  described  completely  in  Appendix  A,  which 
gives  the  flow  diagram,  a  format  guide  for  data  input  cards,  and  the  de¬ 
tailed  program  listing.  As  shown,  both  the  flow  diagram  and  the  program 
listing  contain  some  elements  that  apply  only  to  the  simulations  in  Sec¬ 
tion  2.4.  These  are  indicated  by  broken  lines  in  the  flow  diagram,  and 
by  C*  with  verbal  instructions  in  the  program  listing. 

The  format  guide  should  enable  a  person  who  has  only  slight  acquaint¬ 
ance  with  computers  to  feed  data  into  this  program.  Not  all  of  the  10 
types  of  data  cards  are  used  for  any  one  run,  of  course.  When  the  orig¬ 
inal  data  is  submitted,  card  9  is  included  but  left  blank  (zero  values) . 
For  successive  runs,  card  9  will  contain  KT  (sequence  number)  and  the 
values  of  IR,  IF,  NOP,  and  FOP  supplied  by  the  computer. 

The  form  of  the  program’s  output  is  described  in  the  next  two  sec¬ 
tions,  and  is  shown  in  Appendix  A  as  the  entire  printout  for  the  produc¬ 
tion  example. 
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2.4  COMPUTER  SIMULATION  EXAMPLES 


The  capabilities  of  the  performance  index  method  were  demonstrated 
by  computing  five  simulation  examples  with  the  PIM  program.  The  same 
hypothetical  data  were  used  for  all  of  the  examples,  but  each  displayed  .• 
a  variation  of  procedure  or  of  some  operating  condition. 

The  detailed  description  of  Example  1  indicates  the  general  patterns 
of  input  and  output  for  all  of  the  simulation  cases.  The  shape  of  each 
optimization  analysis  and  its  results  are  shown  graphically  in  Figs.  4 
to  9. 

ASSUMED  DATA 

Operation  and  Tooling.  The  hypothetical  job,  tooling,  and  material 
for  the  simulation  examples  were  taken  to  be 

Operation  Straight  turning 

Machine  tool  Engine  lathe  (IDENT  =  1) 

Workpiece  Custom-455  steel,  6  in.  diam,  24  in.  long 

Performance  Index.  The  same  performance  index  was  used  for  all  of 
the  analyses  —  cost  per  piece,  Q  =  0.  To  evaluate  this  index  by  means 
of  the  subroutine  described  under  COMPUTER  SIMULATION  TESTS  in  Section 
2.2,  the  following  data  were  assumed: 

Time  data  Tool  changing  time 

Machine  handling  time 

Cost  data  Labor  and  overhead  rate 

Cutting  tool  cost 

Machine  data  Maximum  horsepower 

Specific  horsepower 
Machine  efficiency 

Tool  life  data  VT°  *  2  0  3FQ  ' 1 9 *t  =  207 

Constraints .  In  addition  to  the  constraints  imposed  by  horsepower 
and  machine  efficiency,  the  machine  tool  had  assumed  limits  of  capability: 

Available  spindle  speed  range  20  to  1000  rpm  in  50  steps 
Available  feed  range  0.0011  to  0.0168  ipr  in  24  steps 

The  given  machining  operation  imposed  further  restrictions: 

Usable  speed  range  200  to  1000  fpm  in  26  steps  (127  to  637  rpm) 
Usable  feed  range  0.0051  to  0.0102  ipr  in  9  steps 

SIMULATION  EXAMPLE  1 

This  is  a  simple  case  in  which  the  optimum  condition  is  found  by  only 
two  steps  of  optimization  analysis.  As  indicated  by  the  other  examples, 
more  steps  are  usually  required.  The  analysis  of  this  example  is  illus¬ 
trated  in  Fig.  4. 


1  min 
5  min 

$0 . 10/min 
$0. 20/edge 

7.5 

0.75/in. 3/min 
60% 
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The  variants  for  Example  1  were  as  follows: 

Initial  test  points  Selected  by  computer 

Feed  range  Roughing  condition,  0.0051  to  0.0102  ipr 

Error  added  5% 

In  all  of  these  simulations,  a  cost  per  piece  function  was  used  to 
represent  actual  production  costs.  To  better  approximate  shop  results, 
an  arbitrary  error  percentage,  usually  5%  but  in  one  case  10%,  was  added 
to  the  computed  costs.  That  is,  the  indicated  percentage  of  the  cost  was 
taken  and  multiplied  by  a  random  number  that  was  generated  from  -1  to  +1, 
and  this  random  error  was  then  added  to  or  subtracted  from  the  cost  cal¬ 
culated  by  the  cost  per  piece  function. 

Data  Input.  Since  the  computer  is  going  to  determine  the  initial 
tests  points,  we  need  the  following  data  input  to  the  program.  The  code 
terms  are  as  defined  in  Section  2.3,  and  the  numerical  data  are  from  our 
general  assumptions  for  the  simulations  and  the  variant  conditions  given 
for  Example  1. 


IDENT 

= 

1 

M 

24 

VLOW 

- 

127.0 

F  (I)  = 

0.0011  ' 

VHIGH 

= 

637.0 

to 

0.0168 

24  steps 

FLOW 

— 

0.0051 

KT 

■N 

0 

FHIGH 

- 

0.0102 

IR 

0 

VFLOW 

200.0 

One  blank  card 

NPT  = 

0 

for  these  five 

VFHIGH 

= 

1000.0 

VMAX  - 

0.0 

zero  values 

SHP 

= 

0.75 

FMAX  = 

0.0 

HPM 

= 

7.5 

VL(I) 

EFF 

=■. 

0.60 

FD(I) 

T  (I) 

DC 

= 

0.1 

Required  only  when 
NPT  ^  0,  as  in 

RLO 

= 

0.1  [C  ] 
o 

P(D 

Example  3 

TLC 

= 

0.2  [C  ] 
e 

TP  (I) 

Q 

= 

0. 

D 

= 

6 . 

VI 

= 

20.0 

VA 

= 

1000.0 

VCH 

=5 

20.0 

The  information  for  evaluating  the  performance  index  in  the  simula¬ 
tion  program  is  supplied  by  a  data  statement  in  the  subroutine  MACSIM. 
The  statement  reads: 


DATA' A, BB, C,TCT, CO, CE,TM,D, XL/0. 203, 0.194,270.,!. ,0.1, 0.2, 5. ,6. ,24./ 
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FIG.  4.  Computer  simulation  Example  1,  optimization 
analysis.  Cost  per  piece  index  with  5%  error,  initial 
test  points  selected  by  computer,  roughing  condition. 


in  which 


A  =  alpha,  exponent  of  tool  life  =  0.203 

BB  =  beta,  exponent  of  feed  =  0.194 

C  =  constant  in  tool  life  equation  =  270. 

TCT  =  tool  changing  time  =  1. 

CO  =  labor  and  overhead  rate  =  0.1 

CE  =  tool  cost  per  edge  =0.2 

TM  =  machining  time  =  5. 

D  =  diameter  of  workpiece  =  6. 

XL  =  length  of  workpiece  =  24. 


The  error  percentage  is  set  by  a  card  in  the  subroutine  as:  ERR  =  0,05. 

Initial  Test  Points.  The  computer  selected  six  initial  test  points 
for  Example  1: 


itting  Speed 

Feed , 

rPm 

fpm 

ipr 

240 

377 

0.0068 

240 

377 

0.0102 

240 

377 

0.0092 

380 

597 

0.0068 

380 

597 

0.0078 

460 

723 

0.0068 
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The  location  of  these  points  is  shown  by  the  open  triangles  in  Fig.  4. 
The  program  is  expected  to  choose  points  in  a  geometric  arrangement,  but 
in  this  case  the  choice  was  influenced  by  the  horsepower  constraint. 

First  Optimization.  The  printout  for  the  first  optimization  of  the 
simulated  data  read: 


Cutting  Speed 

Feed, 

Unit  Cost 

Prod.  Rate 

rpm 

fpm 

iPr 

Index,  CU 

Index,  PR 

240 

377 

0.0068 

2.169 

0.05 

240 

377 

0.0092 

1.813 

0.06 

240 

377 

0.0102 

1.697 

0.07 

380 

597 

0.0068 

2.611 

0.06 

380 

597 

0.0078 

2.371 

0.06 

460 

723 

0.0068 

3.887 

0.05 

Based  on  these 

data, 

the  program  recommended  the  following 

condition  as  the  optimum  conditions  to  be  tested  further: 


Cutting  speed  =  260  rpm  or  408.4  fpm 
Feed  =  0.0102  ipr 

Further  Optimization.  To  refine  the  results  obtained  with  the  ini¬ 
tial  test  points,  the  computer  then  recommended  four  more  points  for 
testing : 


itting  Speed 

Feed , 

rpm 

_iE L 

140 

220 

0.0102 

260 

408 

0.0102 

380 

597 

0.0078 

500 

785 

0.0060 

The  simulation  program’s  analysis  of  these  points  gave  the  following  data 


Cutting  Speed 
rpm  fpm 


Feed , 

--iP-F- 


Unit  Cost  Prod.  Rate 

Index,  CU  Index,  PR 


140 

220 

0.0102 

2.223 

0.04 

260 

408 

0.0102 

1.599 

0.07 

380 

597 

0.0078 

2.509 

0.06 

500 

785 

0.0060 

4.718 

0.04 

For  these  data,  the  recommended  optimum  conditions  were  again 


Cutting  speed  =  260  rpm  or  408.4  fpm 
Feed  =  0.0102  ipr 

Stopping  Criterion.  Since  the  optimum  conditions  proved  to  be  the 
same  for  the  second  run  as  for  the  first,  the  stopping  criterion  had  been 
met  and  no  further  testing  was  indicated.  For  the  given  machining  opera¬ 
tion,  the  heavily  circled  point  in  Fig.  4  is  the  final  optimum. 
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SIMULATION  EXAMPLE  2 


The  only  difference  between  this  example  and  Example  1  is  that  the 
error  level  was  raised  to  10%.  Because  the  increased  error  caused  more 
oscillation  of  the  optimum  points,  finding  the  final  optimum  required 
three  steps  of  analysis.  Nevertheless,  the  result  for  Example  2,  shown 
in  Fig.  5,  was  the  same  as  that  for  Example  1. 

SIMULATION  EXAMPLE  3 

In  all  other  respects  the  same  as  Example  1,  this  simulation  demon¬ 
strates  selection  of  the  initial  test  points  by  the  user. 

For  Example  3,  nine  initial  points  were  chosen  at  an  arbitrary  loca¬ 
tion  within  the  given  limits  of  the  constraints.  To  rule  out  computer 
selection,  the  following  data  were  added  to  the  original  input : 

NPT  =  9 

VL(I)  =  3  values 

FD(I)  =  3  values  ^  cards 

The  cutting  speeds  chosen  for  this  simulation  were  140,  150,  and  160  fpm. 
The  feeds  were  0.0051,  0.0060,  and  0.0078  ipr .  For  each  VL(I)  there  were 
three  cards  combining  that  value  with  each  of  the  three  FD(I)  values. 

With  these  nine  cards  in  the  deck,  the  computer  program  was  ready  to  be¬ 
gin  the  optimization  procedure. 

Because  of  some  degree  of  oscillation  of  the  optimum  points  in  the 
intermediate  steps,  reaching  the  final  optimum  took  five  steps  of  opti¬ 
mization  analysis.  As  shown  in  Fig.  6,  the  PIM  program  can  eventually 
find  the  optimum  point  regardless  of  how  the  initial  test  points  were 
chosen.  But  manual  selection  is  likely  to  take  more  test  steps. 

SIMULATION  EXAMPLE  4 

To  test  the  versatility  of  the  PIM  program,  in  Example  4  the  initial 
test  points,  again  manually  chosen,  were  located  far  from  the  expected 
optimum  points,  and  the  feed  range  was  broadened  to  0.003  to  0.0102  ipr, 
which  would  cover  both  roughing  and  finishing  conditions. 

Figure  7  shows  that  the  program's  analysis  of  this  simulation  arrived 
at  the  same  recommended  optimum  conditions  as  in  all  previous  cases,  260 
rpm  or  408.4  fpm  and  0.0102  ipr.  It  shows  also  that  five  steps  of  anal¬ 
ysis  were  needed  to  get  there.  Although  the  optimum  point  can  be  found 
with  any  set  of  initial  test  points,  the  closer  the  test  points  are  to 
the  expected  optimum  points,  the  shorter  the  analysis. 

SIMULATION  EXAMPLE  5 

For  this  example  we  return  to  nine  initial  test  points  selected  by 
the  computer  program.  The  purpose  of  Example  5  is  to  demonstrate  how  the 
program  works  with  a  different  machining  constraint,  the  feed  range 
reduced  to  0.003  to  0.0051  ipr  for  finishing  conditions. 
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FIG.  5.  Computer  simulation  Example  2,  optimization 
analysis.  Cost  per  piece  index  with  10%  error,  initial 
test  points  selected  by  computer,  roughing  condition. 


CUTTING  $P££D ,  FPm 


FIG.  6.  Computer  simulation  Example  3,  optimization 
analysis.  Cost  per  piece  index  with  5%  error,  manual 
selection  of  initial  test  points,  roughing  condition. 
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FIG.  7.  Computer  simulation  Example  4,  optimization 
analysis.  Cost  per  piece  index  with  5%  error,  manual 
selection  of  initial  test  points  far  from  expected 
optimum,  feed  rate  extended  to  include  finishing. 


CUTTING  SPCEO,  PPM 


FIG.  8.  Computer  simulation  Example  5,  optimization 
analysis.  Cost  per  piece  index  with  5%  error,  initial 
test  points  selected  by  computer,  feed  range  restrict¬ 
ed  to  finishing  condition. 
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The  results  for  this  simulation  are  shown  in  Fig.  8.  Note  that  be¬ 
cause  of  the  lower  feed  range  the  computer  did  not  have  to  dodge  the 
horsepower  constraint  as  it  did  in  Examples  1  and  2.  It  chose  the  nine 
initial  points  in  the  desirable  geometric  grid.  This  case  took  four  steps 
of  optimization  analysis  to  reach  the  final  recommendation  of  300  rpm  or 
471.2  fpm  and  0.0051  ipr . 


2.5  PRODUCTION  APPLICATION 


When  the  PIM  program  is  used  to  optimize  an  actual  production  case, 
all  of  the  initial  input  is  real  shop  data  and  each  test  point  proposed 
by  the  computer  is  actually  tested  in  the  shop.  The  results  of  the  ma¬ 
chining  tests,  fed  back  to  the  program,  are  the  basis  for  the  program's 
stepwise  estimates  of  optimums  and  its  selection  of  additional  points  to 
be  tested  until  the  final  optimum  point  is  found. 

Described  in  this  section  is  a  "production  run"  of  the  PIM  program, 
with  the  Penn  State  Machinability  Laboratory  and  its  facilities  serving 
as  the  shop. 

NOTES  ON  PROCEDURE 

Program  Adjustments.  As  the  PIM  program  is  shown  in  Appendix  A,  it 
includes  some  cards  that,  apply  only  to  the  simulation  procedure  for  Sec¬ 
tion  2.4.  For  a  production  run,  all  cards  with  C*  in  the  first  two  col¬ 
umns  must  be  converted  to  active  FORTRAN  IV  statements <  That  is  done  by 
reproducing  those  cards,  leaving  blanks  in  place  of  C*.  Also,  any  card 
marked  for  removal  from  the  deck  must  be  removed.  When  these  two  adjust¬ 
ments  are  made,  the  card  deck  will  be  ready  to  run  with  real  shop  data. 

Since  cost  data  are  supplied  in  the  initial  input  and  tool  life  and 
time  data  are  generated  by  the  shop  tests,  the  MACSIM  subroutine  used  in 
the  simulations  is  not  needed  in  production  applications. 

Selection  of  Initial  Test  Points.  One  of  the  basic  functions  of  the 
PIM  program  is  to  make  an  efficient  choice  of  initial  test  points.  It 
performs  this  function  in  the  production  example. 

The  user  may,  however,  choose  the  initial  test  points  himself  (Sec¬ 
tion  2.4,  Example  3),  In  that  case  he  must  collect  all  of  the  information 
needed  for  determining  the  proper  test  values.  He  must  know  what  con¬ 
straints  will  affect  the  choice,  because  he  cannot  pick  points  outside 
these  constraints.  He  must  know  what  spaed  and  feed  settings  are  availa¬ 
ble.  And  he  must  be  sure  to  choose  at  least  six  points,  preferably  as 
many  as  nine,  spaced  several  discrete  settings  apart  and  arranged  in  a 
geometric  pattern  (Section  2.2,  Step  1). 

If  the  machine  tool  to  be  used  is  stepless,  artificial  steps  will 
have  to  be  calculated  because  PIM  computation  depends  upon  discrete  values 
of  the  speed  and  feed  parameters.  The  computer  program  will  calculate 
these  steps  if  the  speed  range  and  the  desired  step  size  are  designated 
in  the  data  deck.  But  the  user,  if  he  wishes,  may  calculate  the  steps 
and  read  them  into  the  deck. 


26 


PRODUCTION  TESTS 


However  the  initial  test  points  are  chosen,  the  result  will  be  six 
or  more  sets  of  cutting  speed  and  feed  settings,  each  of  which  must  be 
tested  on  the  machine  tool.  As  the  optimization  analysis  proceeds,  the 
computer  will  recommend  additional  points  to  be  tested  in  the  shop.  To 
avoid  the  effects  of  error,  which  show  as  oscillations  of  the  optimum, 
it  is  important  to  keep  the  test  conditions  as  consistent  as  possible. 

Each  test  should  be  run  for  the  longest  period  of  time  that  is  prac¬ 
tical.  The  time  period  of  a  test  should  be  at  least  several  times  greateir 
than  both  the  tooling  time  for  each  piece  and  the  life  of  the  cutting 
tool.  The  periods  do  not  have  to  be  equal,  but  no  one  test  period  should 
be  conspicuously  shorter  than  the  others.  If  some  unusual  event,  such  as 
a  long  machine  breakdown  or  an  operator  delay,  interrupts  the  test  period, 
the  lost  time  should  be  subtracted  from  the  total  time  for  the  test. 

Brief  interruptions  that  are  normal  in  production  operations  should  not 
be  deducted. 

For  each  set  of  speed  and  feed  settings  tested,  the  operator  will 
return  three  numbers  to  the  program  for  further  analysis: 

T(I)  =  time  period  of  test,  min 
P(I)  =  number  of  pieces  produced 
TP (I)  =  number  of  tool  changes 

When  a  piece  is  not  completed  or  a  cutting  tool  is  still  in  use  at  the 
beginning  or  end  of  the  period,  the  piece  or  the  tool  can  be  either 
counted  or  left  out,  according  to  the  operator's  judgment. 

PRODUCTION  EXAMPLE:  TEST  CONDITIONS 

To  put  the  PIM  program  to  a  realistic  production  test,  the  initial 
input  information  was  taken  from  our  "shop,"  the  Machinability  Laboratory. 
Computer  runs  were  made  at  each  step  just  as  they  would  be  in  a  plant  ma¬ 
chine  shop,  and  the  tool  life  tests  were  performed  on  a  laboratory  machin¬ 
ability  lathe. 

Operation  and  Tooling.  The  job,  tooling,  and  material  for  the  pro¬ 
duction  example  were  as  follows: 

Operation  Straight  turning 

Machine  tool  Engine  lathe  (IDENT  =  1) 

Workpiece  Custom-455  steel,  6  in.  diam,  24  in,  long 

Cutting  tool  C-2  carbide  insert 

Performance  Index.  As  for  the  simulation  examples,  the  performance 
index  chosen  was  cost  per  piece,  Q  =  0.  The  cost  data  determined  for  this 
particular  shop  were 

Labor  and  overhead  rate  $0.10/min 

Cutting  tool  cost  $0. 20/edge 
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Machining  Constraints.  These  constraints  were  determined  by  the 
characteristics  of  the  lathe: 

Maximum  horsepower  7.5 
Specific  horsepower  0 . 75/in . 3 /min 
Machine  efficiency  60% 

(Note  that  the  horsepower  constraint  is  not  essential  to  the  program.  It 
can  be  optional  in  plant  applications,  depending  on  the  availability  of 
the  data.) 

Available  spindle  speed  range  20  to  1000  rpm  in  50  steps 
Available  feed  range  0.0011  to  0,0168  ipr  in  24  steps 

From  past  experience,  it  was  possible  to  estimate  for  this  material 
and  this  carbide  cutting  tool  the  practical  operating  range  for  rough  ma¬ 
chining  : 

Usable  speed  range  200  to  500  fpm  (127  to  320  rpm) 

Usable  feed  range  0.0051  to  0.0102  ipr 
Depth  of  cut  0.10  in. 

PRODUCTION  EXAMPLE:  INITIAL  DATA  INPUT 


Having  gathered  these  data,  we  are  now  ready  to  instruct  the  computer 
program  to  select  the  initial  test  points.  The  following  data  are  read 
into  the  card  deck  (see  Section  2.3  for  code  definitions): 


o 


IDENT 

= 

1 

VLOW 

= 

127.0 

VHIGH 

as 

320.0 

FLOW 

= 

0.0051 

FHIGH 

= 

0.0102 

VFLOW 

= 

200.0 

VFHIGH 

= 

500.0 

SHP 

= 

0.75 

HPM 

= 

7.5 

EFF 

= 

0.60 

DC 

= 

0.1 

RLO 

= 

0.1 

TLC 

0.2 

Q 

= 

0. 

D 

= 

6. 

VI 

= 

20.0 

VA 

as 

1000.0 

VCH 

— 

20.0 

M 

F  (I) 


KT 

IR 

NPT 

VMAX 

FMAX 


=  24. 


=  0.0011  ] 


to 


24  steps 


0.0168  J 


=  0 
=  0 
=  0 
=  0 
=  0 


Zero  on  initial  run. 
Computer  supplies  values 
for  succeeding  runs. 


VL(I)  ' 
FD(I) 
T(I) 
P(I) 

TP (I)  , 


No  cards  needed  for 
initial  run.  NPT  cards 
needed  for  succeeding 
runs . 
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PRODUCTION  EXAMPLE:  PIM  ANALYSIS 

The  complete  computer  printout  for  the  production  example  is  given 
in  Appendix  A  and  the  analysis  is  shown  graphically  in  Fig.  9. 

Because  this  example  is  intended  to  demonstrate  how  the  program  oper¬ 
ates  in  a  production  application,  the  analytic  process  is  described  as  a 
stepwise  procedure  with  some  general  instructions  and  notes.  The  reader 
will  find  it  useful  to  follow  these  steps  in  the  computer  printout,  remem¬ 
bering  that  all  numerical  data  apply  to  the  example  only. 

Step  1:  Initial  Test  Points  and  Tests.  Based  on  the  data  input,  the 
first  output  of  the  program  is  the  nine  initial  test  points,  shown  in  the 
three  columns  headed  "Test  Conditions"  in  the  following  tabulation.  These 
are  actually  combinations  of  three  speeds  and  three  feeds.  For  a  time 
period  of  240  min  (4  hr)  in  this  case,  each  combination  was  tested  on  the 
lathe,  with  these  results: 


Test 

Conditions 

Production  Test 

Results 

Cutting  Speed 

Feed, 

Time , 

No .  of 

No.  of 

rpm 

fpm 

ipr 

min 

Pieces 

Tools 

160 

251 

0.0068 

240 

11 

1 

160 

251 

0.0092 

240 

15 

2 

160 

251 

0.0102 

240 

16 

2 

220 

346 

0.0068 

240 

15 

7 

220 

346 

0.0092 

240 

20 

9 

220 

346 

0.0102 

240 

22 

10 

260 

408 

0.0068 

240 

18 

16 

260 

408 

0.0092 

240 

24 

21 

260 

408 

0.0102 

240 

27 

23 

The  results  are  fed  back  to  the  computer  program  by  adding  the  fol¬ 
lowing  cards : 

NPT  =  9  (Number  of,  tests  points  run,  IF) 


VL  ( I. ) 
FD(I) 
T  (I) 
P(I) 
TP  (I) 


Settings  and  results 
for  each  test  (9  cards) 


Note  that  in  real  production  the  operator  might  have  chosen,  from  his 
own  experience,  to  eliminate  the  group  of  tests  at  the  lowest  cutting 
speed  (160  rpm) .  He  would  then  report  only  the  remaining  six  tests,  as 
NPT  =  6.  But  at  least  six  tests  should  be  fed  back  for  the  regression 
analysis. 


Note  also  that  the  test  runs,  though  they  vary  in  efficiency,  are  not 
"lost  time"  in  production.  During  the  36  hr  spent  on  the  first  tests  for 
this  example,  168  pieces  were  actually  produced. 
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FIG.  9.  Production  example  with  shop  test  data. 
Cost  per  piece  index,  initial  test  points  select 
ed  by  computer,  roughing  condition. 


Step  2:  First  Optimization.  At  the  end  of  the  first  run,  the  com¬ 
puter  prints  two  codes  for  the  next  run.  For  this  example  they  are 


IR  =  0 
IF  =  9 

The  IR  code  indicates  whether  the  test  results  will  be  fitted  to  a  sur¬ 
face  regression  analysis  (IR  =  0)  or  to  a  curvilinear  regression  analysis 
(IR  =  1) .  The  IF  code  indicates  the  number  of  test  results  (NPT)  to  be 
reported.  Along  with  the  Step  1  test  data,  these  codes  must  be  returned 
to  the  computer  before  the  first  optimization  run  of  the  program. 

At  this  point  we  also  activate  the  KT  code  to  indicate  the  sequence 
number  of  this  analysis,  as  KT  =  1.  By  raising  the  KT  number  with  each 
succeeding  run  (KT  =  2  and  so  on)  we  get  numbered  printouts  that  can 
prevent  confusion,  particularly  when  the  number  of  optimization  runs  is 
greater  than  two.  Use  of  the  KT  numbers  is  optional,  however;  the  pro¬ 
gram  can  be  run  without  them. 

The  second  computer  run  produces  the  first  optimum  condition  for  the 
chosen  performance  index,  cost  per  piece: 

Cutting  speed  320  rpm  or  502.7  fpm 
Feed  0.0094  ipr 
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It  then  recommends,  in  this  case,  three  conditions  for  further  test¬ 
ing  (shown  in  the  following  tabulation).  One  of  these  conditions  will 
always  he  the  optimum  condition  previously  computed.  The  printout  again 
ends  with  the  codes  for  the  next  run. 

The  three  new  points  were  tested  on  the  lathe,  with  due  concern  for 
the  cautions  indicated  in  the  paragraphs  headed  PRODUCTION  TESTS.  These 
tests  yielded  the  following  feedback  to  the  computer  (three  cards) : 


Test  Conditions  Production  Test  Results 


Jtting  Speed 

Feed , 

Time , 

No.  of 

No .  of 

rpm 

fpm 

ipr 

min 

Pieces 

Tools 

240 

377 

0.0094 

240 

23 

14 

280 

440 

0.0094 

240 

26 

31 

320 

503 

0.0094 

240 

30 

59 

Step  3:  Second  Optimization.  At  the  start  of  this  run  (KT  -  2)  we 
must  again  pick  up  the  codes  printed  at  the  end  of  the  Step  2  computer 


output : 


IR  =  1 
IF  =  3 

This  means  that  the  regression  analysis  of  the  three  new  test  points  will 
be  curvilinear. 

On  the  next  computer  printout  we  get  the  second  optimum  point  recom¬ 
mendation: 


Cutting  speed  300  rpm  or  471.2  fpm 
Feed  0.0102  ipr 

This  time  the  computer  recommends  four  conditions  for  additional 
testing,  one  of  which  is  the  second  optimum.  The  codes  given  for  the 
next  run  are  IR  =  1,  IF  =  4. 

The  tests  were  run  on  the  lathe,  and  the  following  results  were  re¬ 
ported  to  the  program: 


Test  Conditions 


Cutting  Speed 

Feed, 

rpm 

fpm 

lP.r— 

220 

346 

0.0102 

260 

408 

0.0102 

300 

471 

0.0102 

320 

503 

0.0094 

Note  that  the  time  period  of 

:  these 

consistent . 


Production  Test 

Re  suits 

Time , 

No .  of 

No .  of 

min 

Pieces 

Tools 

246 

22 

10 

246 

27 

23 

246 

31 

47 

246 

30 

59 

was  increased,  but  it  was  kept 


Step  4:  Third  Optimization.  Carding  in  the  codes  from  the  end  of 
the  previous  printout  and  the  new  data  from  Step  3,  we  again  run  the 
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program.  The  result  is  the  third  recommended  optimum: 

Cutting  speed  300  rpm  or  471  fpm 
Feed  0.0102  ipr 

Stopping  Criterion.  The  optimum  conditions  computed  by  the  last  two 
runs  are  identical.  The  stopping  criterion  (Section  2.2,  Step  6)  has 
therefore  been  satisfied,  and  no  further  testing  is  indicated.  If  the 
stopping  criterion  had  not  yet  been  reached,  the  computer  would  have 
recommended  a  new  set  of  test  points  and  the  procedure  would  continue  as 
before . 

It  is  important  to  note  that  the  program  user  may  exercise  certain 
options  as  the  computations  approach  the  final  optimum.  If  the  last  two 
optimums  are  only  one  or  two  settings  apart,  he  may  choose  to  end  the 
procedure  before  the  stopping  criterion  is  met.  The  true  optimum  is  prob¬ 
ably  on  the  continuum  between  those  settings,  and  any  of  the  two  or  three 
settings  may  be  adequate  for  the  accuracy  desired.  If  the  estimated  opti¬ 
mum  point  oscillates  among  several  close  values,  it  may  be  reasonable  to 
decide  that  any  of  these  values  is  usable  and  the  testing  can  be  stopped. 

CONCLUSION 

We  chose  for  this  example  a  case  for  which  the  optimum  cutting  con¬ 
ditions  had  been  established.  The  final  optimum  computed  by  the  PIM  pro¬ 
gram  is  very  close  to  the  known  speed  and  feed  figures,  certainly  close 
enough  to  support  confidence  that  the  program  can  be  applied  successfully 
to  actual  production  if  the  test  conditions  are  properly  controlled. 
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PART  3-  PRODUCTION  OPTIMIZATION  METHOD 


When  tool  life  data,  cost  data,  and  time  study  data  are  available, 
the  production  optimization  method  (POM)  provides  a  simple  and  practical 
procedure  for  optimizing  machining  parameters.  POM  requires  fewer  shop 
tests  than  the  performance  index  method  (PIM) ,  and  it  yields  both  more 
accurate  optimization  and  a  larger  quantity  of  useful  information.  It  is 
less  versatile  than  PIM  in  that  it  is  designed  to  analyze  only  turning 
operations,  but  it  offers  the  industrial  user  greater  flexibility  in  the 
selection  of  optimum  conditions  for  metal  cutting.  For  these  reasons, 

POM  is  recommended  as  the  better  method  for  practical  applications  where 
the  necessary  data  are  obtainable. 

Like  the  PIM  program,  the  POM  computer  program  is  constructed  to 
serve  the  machinist  on  the  job.  The  mathematical  processes  described  in 
Section  3.2  are  written  into  it  and  will  operate  automatically  unless  they 
are  modified  by  engineering  decision.  Output  from  the  program  is  in  the 
form  of  printed  tables  giving  a  ranked  range  of  cutting  conditions  for 
roughing  and  finishing,  keyed  to  specific  production  objectives. 

Section  3.3  describes  the  program  structure  and  the  language  code  for 
input  data.  A  sample  application  of  POM  is  presented  in  Section  3.4. 


3.1  CONCEPT  AND  TECHNIQUES 


The  production  optimization  method  analyzes  optimum  cutting  condi¬ 
tions  by  gradual  refinement  of  tool  life  data  fed  back  to  the  program  from 
production  tests  in  the  workshop.  By  an  exhaustive  enumeration  procedure 
the  program  computes  and  ranks  all  combinations  of  cutting  speed,  feed, 
and  depth  of  cut  within  the  constraints  imposed  by  a  given  workpiece,  ma¬ 
chine  tool,  cutting  tool,  and  surface  finish  requirement.  It  also  orders 
these  values  in  terms  of  one  or  more  of  three  objective  functions:  mini¬ 
mum  cost  per  piece,  maximum  production  rate,  and  maximum  profit  rate  (4). 

TOOL  LIFE  DATA 


The  POM  program  has  the  unique  characteristic  of  generating  the  re¬ 
quired  tool  life  data  for  any  given  material  and  cutting  tool.  This  capa¬ 
bility  leads  to  a  high  degree  of  accuracy  and  flexibility  in  production 
optimization.  It  is  also  a  distinct  advantage  in  analyzing  new  workpiece 
and  cutting  tool  materials. 


POM  employs  the  conventional  concept  of  Taylor's  tool  life  equation, 
in  the  general  form  of 


VTaFeDY  = 


C 


The  analysis  begins  with  an  educated  initial  estimate  of  the  tool  life 
parameters  for  the  specific  case,  based  upon  published  tool  life  data  or 
practical  experience,  or  both.  The  values  computed  with  this  estimate 
are  tested  in  production,  and  the  test  results  are  fed  back  to  the  program 
as  data  for  further  computation  to  refine  the  parameter  values.  This 
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process  is  repeated  until  the  tool  life  equations  resulting  from  two  suc¬ 
cessive  iterations  show  no  difference,  or  a  percentage  of  difference  that 
falls  within  a  specified  limit.  The  accepted  tool  life  equation  then 
functions  in  the  POM  program’s  optimization  procedure  for  the  case  being 
analyzed. 

The  production  example  in  Section  3.4  demonstrates  that  inaccuracies 
in  the  estimated  tool  life  parameters  are  cleared  as  the  analysis  prog¬ 
resses.  If  the  initial  guess  proves  to  be  exact  or  if  the  true  tool  life 
parameters  are  known  at  the  outset  and  read  into  the  equation,  the  program 
will  quickly  verify  the  values  and  go  to  the  optimization  procedure. 

The  mathematics  involved  in  computing  the  tool  life  equations  is  ex¬ 
plained  in  Section  3.2.  The  intention  at  this  point  is  only  to  indicate 
the  efficiency  and  simplicity  of  POM  treatment  as  compared  with  the  usual 
means  of  obtaining  tool  life  information.  Because  every  production  ma¬ 
chining  job  is  influenced  by  plant  factors,  published  tool  life  data, 
though  abundant,  can  provide  no  more  than  approximations  for  specific 
cases.  For  new  workpiece  materials  ®r  new  cutting  tool  materials,  the 
published  data  are  not  likely  to  be  helpful  except  for  rough  estimating. 

When  a  tool  life  equation  must  be  evaluated,  the  usual  resort  is  to 
elaborate,  time-consuming,  and  expensive  tool  life  tests  at  a  machin- 
ability  laboratory.  Such  tests  are  performed  under  closely  controlled 
conditions,  with  tool  life  criteria  based  upon  a  predetermined  limit  of 
either  flank  wear  or  crater  wear.  There  is  no- guarantee  that  the  result¬ 
ing  tool  life  data  can  be  readily  applied  to  production  operations  with  a 
particular  machine  tool  in  a  particular  workshop.  Unless  they  are  well 
tested  and  adjusted  to  accommodate  local  restrictions,  tool  life  data  ob¬ 
tained  in  a  laboratory  often  prove  to  be  inaccurate  in  the  shop  (6, 7,8,9). 

Tool  life  data  computed  by  the  POM  program  are  fitted  to  the  indi¬ 
vidual  case.  The  shop  tests  are  performed  in  the  course  of  normal  pro¬ 
duction,  the  computations  are  automatic  and  rapid,  and  the  results  are 
ready  for  immediate  use. 

MATHEMATICAL  TECHNIQUES 

The  mathematical  technique  for  evaluating  the  tool  life  equation  is 
to  compute  the  exponential  values  and  the  constant  by  solving  simultane¬ 
ous  equations  with  a  minimum  of  data.  At  least  three  sets  of  data  for 
speed,  feed,  and  tool  life  are  needed  to  compute  a,  $,  and  C.  (Since 
depth  of  cut  is  usually  the  least  critical  of  the  tool  life  parameters, 
for  this  analysis  we  assume  that  y  =  0) .  The  initial  equations  are  solved 
with  estimated  values,  which  are  then  refined  by  additions  of  test  data 
from  the  shop. 

Conventionally,  the  method  of  partial  derivatives  is  applied  to 
optimization  analysis  of  machining  conditions.  That  method  can  optimize 
only  one  parameter  at  a  time,  either  cutting  speed  or  feed,  with  the  other 
parameters  held  constant.  It  does  not  present  a  combination  of  all  pa¬ 
rameters  for  the  optimum  condition.  Also,  the  cutting  speed  for  the  maxi¬ 
mum  profit  rate  cannot  be  determined  explicitly  by  this  method  except  for 
particular  values  of  the  tool  life  exponents. 
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POM  analysis  employs  instead  an  exhaustive  enumeration  procedure, 
subject  to  a  series  of  cutting  speeds,  feeds,  and  depths  of  cut.  By 
multiregression  analysis,  the  program  selects  a  range  of  optimum  cutting 
conditions  ranked  according  to  one  of  three  production  objectives.  As 
the  general  program  is  now  written  (Appendix  B) ,  the  final  optimization 
automatically  prints  rankings  for  minimum  cost  per  piece,  maximum  pro¬ 
duction  rate,  and  maximum  profit  rate,  for  both  roughing  and  finishing 
conditions . 


3.2  MATHEMATICAL  PROCEDURES 

This  section  explains  the  calculations  essential  to  the  POM  program, 
which  include  the  creation  of  artificial  steps  for  a  stepless  lathe,  the 
computation  of  cutting  speed  increments,  evaluation  of  the  tool  life  equa¬ 
tion,  and  the  relations  for  the  three  production  objectives.  Since  the 
necessary  procedures  are  written  into  the  program,  the  production  user 
need  not  concern  himself  with  the  mathematics.  The  information  herein 
will  be  needed,  however,  if  the  general  program  is  revised  for  a  closer 
fit  to  the  requirements  of  an  individual  plant. 

ARTIFICIAL  SPINDLE  SPEED  STEPS 

Like  the  PIM  program,  the  POM  computer  program  is  designed  to  opti¬ 
mize  cutting  conditions  for  either  a  stepped  lathe  (IDENT  =  0)  or  a  step¬ 
less  lathe  (IDENT  =  1),  For  a  stepped  lathe,  the  speed  steps  are  simply 
read  in  as  known  values,.  For  a  stepless  lathe,  the  program  requires  arti¬ 
ficial  steps,  which  can  be  calculated  in  the  following  manner. 

The  minimum  and  maximum  spindle  speeds  for  a  given  lathe  are  always 
known.  The  user  will  specify  the  size  of  increment  desired  between  spin¬ 
dle  speed  steps.  With  these  data,  the  number  of  steps  can  be  calculated 
as 

X  =  [(RPM  -  RPM  )/RPM  ,]  +  1.0  (12) 

a  i  ch 

where  X  =  number  of  increments,  a  decimal  number 
RPM  =  highest  spindle  speed  available,  rpm 

RPM^  =  lowest  spindle  speed  available,  rpm 

RPM  ,  =  size  of  increment 
ch 

But  X  expressed  as  a  decimal  number  has  to  be  changed  to  an  integer  L,  as 
L  =  [X].  So  we  add  1  to  X  to  make  the  integer  L,  the  largest  number  in  X, 
large  enough  to  span  the  range  of  spindle  speeds. 

The  possibility  still  exists,  however,  that  L  will  be  less  than  X. 

In  that  case,  the  computer  program  must  automatically  increase  L  by  1. 

The  program  therefore  computes  RPM^,  the  values  of  the  L  steps,  as 

RPM±  =  RPM^  +  (I  -  l)RPMch  I  =  1,2,...L  (13) 
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The  range  of  I  must  begin  with  1  because  the  computer  cannot  index  from 
zero. 


CUTTING  SPEED  INCREMENTS 

Having  set  the  incremental  values  of  the  spindle  speed,  we  can  now 
calculate  the  cutting  speed  increments  and  the  minimum  and  maximum  cutting 
speeds  possible  for  this  operation.  Because  cutting  speed  is  a  function 
of  the  workpiece  diameter,  we  must  knew  the  initial  and  final  diameters 
for  these  calculations: 

V  .  =  7TD  (RPM  )/12.0  (14) 

min  s  i 

where  V  .  =  cutting  speed,  fpm,  at  initial  diameter 

with  lowest  spindle  speed  (rpm) 

Dg  =  initial  diameter  ef  workpiece,  in. 

RPM^  =  lowest  spindle  speed  available,  rpm 

and 

V  =  7TD _ (RPM  )/12.0  (15) 

max  f  a 

where  V  =  cutting  speed,  fpm,  at  final  diameter 
with  highest  spindle  speed  (rpm) 

Df  =  final  diameter  of  workpiece,  in. 

RPM  =  highest  spindle  speed  available,  rpm 
3 


The  cutting  speed  increments  can  then  be  calculated  with  equations 
in  the  form  of 


V  exp  (log  V  /V  ) 
min  max  min 


(K-l) / (L-l) 


K  =  1,2, ...L 


This  is  a  step  function  used  to  form  L  steps  between  V  and  .  Since 

exp  (log  X)  =  X,  we  can  write  the  equations  as 


V  =  V  .  (V  /V  .  )(K 
k  min  max  min 


K  =  1,2, ...L 


(16) 


It  follows  that  if  K  =  1, 

V  =  V  .  (V  /V  ,  ) 0  =  V  , 
l  min  max  min  min 

and  if  K  =  L, 

V.  =  V  .  (V  /V  ,  )*  =  V 
L  min  max  min  max 
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PARAMETERS  OF  TOOL  LIFE  EQUATION 


Since  depth  of  cut  appears  to  be  the  least  significant  factor  for 
this  analysis,  we  can  reduce  the  Taylor  tool  life  equation  to 

VTaF3  =  c  (17) 

where  V  =  cutting  speed,  fpm 
T  =  tool  life,  min. 

F  =  feed,  ipr 
a,  3  =  exponents 
C  =  constant 

With  this  as  our  basic  relation,  we  begin  the  iterative  process  of  com¬ 
puting  and  refining  the  tool  life  parameters  a,  3,  and  C. 

Four  situations  occur  in  every  search  for  an  optimum  cutting  condi¬ 
tion:  1)  no  shop  data  are  available,  2)  one  set  of  shop  data  is  availa¬ 

ble,  3)  two  sets  of  shop  data  are  available,  and  4)  three  or  more  sets  of 
shop  data  are  available.  The  program  deals  with  these  situations  stepwise 
as  they  occur. 

Step  1:  No  Shop  Data.  To  obtain  an  initial  set  of  cutting  speeds 
and  feeds,  the  program  is  run  without  any  shop  data  (NDATA  =  0).  The 
tool  life  equation  for  this  step  is 

VTaiF^  =  Ci 

for  which  ai ,  3i,  and  C!  are  estimated  values,  an  "educated  guess,"  sup¬ 
plied  by  the  user. 

With  these  data,  the  program  generates  20  ranked  sets  of  speed  and 
feed  combinations  within  the  capability  of  the  machine  tool.  The  computer 
output  is  a  printed  table  optimized  for  the  given  production  objective  by 
exhaustive  enumeration.  One  set  is  selected  from  the  table  and  tested  in 
the  shop  to  determine  a  tool  life  for  that  combination.  The  user  may 
choose  the  first  set  or  may  exercise  his  own  judgment  in  the  selection. 

The  shop  test  data,  which  should  always  be  an  average  of  several  tests, 
are  returned  to  the  program. 

Step  2;  One  Set  of  Shop  Data.  With  one  set  of  "known"  conditions, 
only  one  of  the  three  unknown  tool  life  parameters  can  be  calculated.  We 
choose  to  calculate  first  the  constant  C.  Calling  the  values  from  the 
shop  tests  Vlt  F: ,  and  T1}  we  find  a  new  value  of  C  from  the  relation 

C  =  V  Ta2F32  (18) 

in  which  012=04  anc*  32  =  3^ 

The  tool  life  equation  for  Step  2  becomes 
VTa2F^2  =  c2 
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and  from  this  equation  the  computer  program  generates  another  20  ranked 
sets  of  speeds  and  feeds.  For  one  of  these  sets,  shop  tests  will  return 
to  the  program  the  data  for  V2,  F2 ,  and  T2.  For  mathematical  reasons  re¬ 
lated  to  the  logarithms  used  later,  care  must  be  taken  to  select  a  test 
set  that  will  yield  values  such  that  V2  ^  Vj  and  F2  ±  F^ 

Step  3:  Two  Sets  of  Shop  Data.  There  are  now  two  distinct  sets  of 
data,  VpFpTj  and  V2,F2,T2,  allowing  two  of  the  three  tool  life  param¬ 
eters  to  be  calculated.  We  choose  C  and  B.  The  new  values  C3  and  B3 
must  satisfy  both  sets  of  data.  That  is, 

V  Ta3F^3  =  C  V  Ta3F^3  =  c 

111  3  2  2  2  3 

in  which  a3  =  a2 .  Since  b®th  of  the  left  sides  equal  C3,  they  must  be 
equal : 

V  T°I'3F^3  =  V  Ta3F^3 

1)1  2  2  2 

or 

(F  /F  )3s  =  (V  /V  ) (T  /T  )a3 
12  2  12  1 

We  take  logarithms  of  both  sides  of  that  equation: 

log  (F  /F  )63  =  log  [(V  /V  ) (T  /T  )a3] 

12  2)21 

By  the  rules  of  logarithms,  the  following  changes  can  be  made: 

B  log  (F  /F  )  =  log  (V  /V  )  +  [a  log  (T  /T  )] 

3  1  2  2  1  3  2  1 

-  -  log  (V  /V  )  -  [a  log  (T  /T  )] 

12.  3  12 

We  then  solve  for  B3: 

B  =  [-  log  V  /V  -  (a.  log  T  ,/T  )  ]  / log  F  /F 

3  12  3  12  12 

Having  obtained  a  value  for  S3,  we  can  solve  for  C3: 

c  =  v  t^f3 3 

3  2  2  2 

The  new  tool  life  equation 
VTa3F33  =  C 

3 

is  used  to  generate  20  more  ranked  combinations  of  speeds  and  feeds. 
Again,  the  one  set  selected  to  test  in  the  shop  for  T3  must  be  such  that 
V 3  ^  V2  ^  Vj  and  F3  ^  F2  4-  Fj,  The  data  from  these  tests  are  returned  to 
the  computer  program  as  V3,F3,T3. 


(19) 


(20) 
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Step  4;  Three  or  More  Sets  of  Shop  Data.  With  at  least  three  sets 
of  data,  we  can  calculate  all  three  of  the  tool  life  parameters.  The 
procedure  described  for  Step  4  can  be  used  only  when  there  are  three  or 
more  sets  of  shop  data  (n  =  3,4,5...). 

By  multiple  regression  analysis,  we  will  seek  a  linear  equation  that 
best  describes  the  discrete  data  points  available.  Taking  the  logarithm 
of  both  sides  of  Eq.  17  and  simplifying  it,  we  get  a  linear  equation  suit¬ 
able  for  regression  analysis: 

log  V  +  (a  log  T)  +  (8  log  F)  =  log  C 

in  which  log  V  is  the  dependent  variable,  log  T  and  log  F  are  the  inde¬ 
pendent  variables,  and  a,  8,  and  C  are  constants.  Rewritten  in  standard 
form,  this  equation  becomes 

log  V  =  log  C  -  (a  log  T)  -  (8  log  F) 

If  log  V  is  taken  to  be  a  value  from  the  desired  equation  and  the  right 
side  an  estimated  value  from  the  data,  the  object  of  this  analysis  is  to 
minimize  the  difference  between  the  two: 

log  V  -  [log  C  -  (a  log  T)  -  (8  log  F) ] 

There  are  n  such  equations,  of  course,  and  the  difference  should  be  mini¬ 
mized  for  all  n  points.  Since  the  differences  D  may  be  plus  or  minus,  we 
will  use  the  sum  of  their  squares: 

n  2 

D  =  l  [log  V  -  log  C  +  (a  log  T  )  +  (8  log  F  ) ] 
i=l  1 


If  the  partial  derivatives  of  D  with  respect  to  the  three  constants 
are  taken  and  set  equal  to  zero,  the  resulting  values  of  the  constants 
will  ensure  the  minimum  possible  value  ®f  D.  That  is, 

n 

3D/(3  log  C)  =  -2  l  [log  V  -  log  C  +  (a  log  T  )  +  (8  log  F  )]  =0 
i=l  1 

n 

3D/3a  =2  l  [log  V  -  log  C  +  (a  log  T.)  +  (6  log  F.)]  log  T.  =0 

i=l  1  1  1  1 

n 

3D/38  =2  l  [log  V.  -  log  C  +  (a  log  T.)  +  (8  log  F.)]  log  F.  =  0 

i=l  1  1  1  1 

Rewriting  the  first  of  these  equations  by  dividing  out  the  -2  and  distrib¬ 
uting  the  summation  signs  gives  . 

n  n  n  n 

l  log  V.  -  l  log  C  +  (a  l  log  T.)  +  (8  l  log  F  )  =  0 

i=l  1  i=l  i»l  1  i=l 
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which  can  be  rearranged  with  £  log  C  =  n  log  C,  so  that 

i=l 

n  n  n 

n  log  C  -  (a  l  log  T  )  -  (0  £  log  F  )  =  £  log  V 

i=l  1  i=l  1  i=l  1 


(21) 


Likewise,  the  second  and  third  equations  can  be  written  as 

n  n  2  n 

log  C  l  log  T  -  [a  l  (log  T  )  ]  -  (0  l  log  F  log  T  ) 
i=l  1  i=l  1  i=l 


-  I  vi  l08  Ti 
i=l  1 


(22) 


n 


n 


log  C  l  log  F  -  (a  l  log  T  log  F  )  -  [0  l  (log  F  )  ] 
i=l  1  i=l  1  i=l  X 


=  l  log  V  log  F 
i=l  1 


(23) 


We  now  have  expressions  for  each  of  the  three  unknowns,  Eqs.  21,  22, 
and  23.  These  equations  can  be  solved  in  any  one  of  a  great  number  of 
ways.  The  POM  program  employs  determinants  to  find  the  values  of  log  C, 
a,  and  0.  The  tool  life  constant  C  is  then  found  by 

C  =  exp  (log  C) 
where  log  C  is  what  is  known. 

The  third  set  of  shop  data  provides  the  first  set  of  calculated  tool 
life  parameters:  a4,04,C4.  When  the  fourth  set  of  shop  data  is  analyzed, 
the  new  parameters  a5,05,C5  are  compared  to  the  first  set  to  determine 
the  percentage  of  change  in  the  values: 

a  =  j (as  -  a4)/a4|  x  100 

3  -  | (05  -  34)/B4|  x  100 

0  =  |  (c5  -  C4)/Cj  x  100 

This  procedure  is  continued  until  the  value  for  each  of  the  three  param¬ 

eters  is  equal  to  or  less  than  the  prescribed  stopping  criterion  for  that 
parameter.  The  resulting  tool  life  equation  is. used  to  generate  the  final 
set  of  solutions. 

PRODUCTION  OBJECTIVES 

Though  job  work  with  POM  will  usually  be  done  in  terms  of  a  single 
production  objective,  the  program  as  it  is  now  written  (Appendix  B) 
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computes  all  three  production  objectives  in  the  final  optimization.  The 
following  are  the  relations  for  calculating  the  objectives. 

1)  Minimum  Coat  per  Piece.  This  computation  involves,  of  course,  all 
of  the  cost,  time,  and  machining  data  supplied  to  the  program.  The  rela¬ 
tion  is 


C  =  C  [N(2T,  +  XL/F  )+Tt+NT  +T]  +  (CT  +C)P 
u  o  1  r  L  gc  gc  c  o  ct  e 


u  o  1  r  l  gc  gc  c  o  ct  e 

where  =  total  machining  cost,  $/piece 

C  =  labor  and  overhead  rate,  $/min 
o 

N  =  number  of  cuts 

Tx  =  time  for  motion  at  start  and  end  of  cut,  min 
XL  =  length  of  workpiece,  in. 

F^  =  return  speed  of  carriage,  in. /min 
T  =  work  handling  time,  min 

Li 

N  =  number  of  gear  changes 
gc 

=  time  per  gear  change,  min 
Tc  =  cutting  time,  min 
T  =  tool  changing  time,  min 
=  tool  cost,  $/edge 
P  =  number  of  edges  per  piece 

2)  Maximum  Production  Rate.  The  total  time  per  piece  Tpro(j,  in  min¬ 
utes,  is  the  sum  of  all  the  time  elements: 

T  ,  =  N(2T,  +  XL/F  )+Tt+NT  +T+TP 
prod  1  r  L  gc  gc  c  ct 

The  production  rate,  expressed  as  60/Tpro(j,  is  the  reciprocal  of  the  time 
per  piece  multiplied  by  60  min,  so  that  it  is  found  as  pieces  per  hour. 

3)  Maximum  Profit  Rate.  This  simple  calculation  relates  all  of  the 
unit  cost  and  time  factors  to  the  price  per  piece: 

Profit  =  (Price  -  C  -  C  )/T 

mt  u  prod 


in  which  Price  is  the  selling  price  per  piece,  C  is  the  material  cost, 
and  the  other  factors  are  as  previously  defined. 


3.3  COMPUTER  PROGRAM  FOR  POM 

For  the  production  optimization  method,  as  for  the  performance  index 
method,  the  computer  program  uses  the  language  of  FORTRAN  IV  and  was  de¬ 
veloped  with  the  IBM  360,  Model  67  computer  in  the  Computation  Center  at 
Penn  State.  The  POM  program  is  adaptable,  however,  to  other  computer 
systems . 
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Since  the  mathematical  operations  explained  in  Section  3.2  are  writ¬ 
ten  into  the  program  and  the  code  terms  for  data  input  are  directly  re¬ 
lated  to  the  vocabulary  of  production,  POM  analysis  can  be  performed  by 
shop  personnel  who  are  not  familiar  with  the  mathematical  techniques  or 
with  computer  programming. 

The  structure  of  the  program  is  described  in  detail  in  Appendix  B, 
which  provides  the  flow  diagram,  the  format  guide  for  data  input  cards, 
and  a  complete  listing  of  the  computer  program.  How  POM  operates  is 
briefly  surveyed  in  this  section  and  demonstrated  by  the  production  exam¬ 
ple  in  Section  3.4. 

For  ready  reference  in  relation  to  the  production  example,  the  input 
data  code  is  listed  and  defined  at  the  end  of  this  section.  The  code 
terms  are  given  in  the  order  of  their  occurrence  in  the  format  guide  in 
Appendix  B. 

GENERAL. PROCEDURE 

The  POM  program  uses  cost  data,  time  study  data,  machine  parameters, 
and  tool  life  data  obtained  in  shop  tests  to  arrive  at  an  optimum  set  of 
cutting  conditions  for  any  one  or  all  of  three  production  objectives.  The 
initial  data,  including  an  estimate  of  tool  life  parameters  for  the  job, 
are  gathered  and  prepared  for  the  computer  according  to  the  format  guide 
(Appendix  B) .  From  this  input  the  computer  generates  with  subroutine 
RANK,  explained  later  in  this  section,  a  ranked  set  of  speeds  and  feeds. 

The  user  selects  one  of  these  combinations  of  feed  and  speed  to  test 
for  tool  .life  in  his  shop.  Note  that  these  tests  are  run  as  normal  pro¬ 
duction  operations,  the  only  requirement  being  that  no  factor  can  be 
changed  during  several  repetitions  of  the  test  for  any  one  combination. 

The  object  is  to  obtain  good  average  figures  for  the  tool  life  parameters, 
which  will  be  fed  back  to  the  program  f«r  refinement  of  the  initial  tool 
life  estimate.  The  results  of  this  shop  test  are  punched  into  a  data  card 
(type  10),  and  this  card  is  stored  behind  the  original  data,  to  be  used  in 
the  second  analytic  run  of  the  program. 

That  procedure  is  repeated,  adding  one  data  card  per  investigation, 
until  the  change  in  each  of  the  tool  life  equation  parameters  is  smaller 
than  some  acceptable  value.  The  final  optimization  analysis  is  then  gen¬ 
erated  with  subroutine  OPTM,  described  later  herein. 

SPINDLE  SPEED  STEPS 

When  all  of  the  job  constants  (cards  1  through  3)  have  been  read  into 
the  program,  the  computer  tests  IDENT  to  determine  the  type  of  spindle 
speed  control  this  machine  tool  has.  If  IDENT  =  0,  the  lathe  is  stepped 
and  the  computer  looks  for  two  or  more  cards,  one  giving  the  number  of 
speed  steps  and  one  or  more  giving  the  values  of  the1 spindle  speed  at  each 
step. 

If  IDENT  =  1,  the  lathe  is  stepless  and  the  computer  looks  for  a 
single  data  card  containing  the  lowest  and  highest  spindle  speeds  availa¬ 
ble  and  the  size  of  the  increment  from  one  speed  to  the  next  (card  4) ,  as 
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specified  by  the  user.  With  these  data,  the  computer  will  calculate  the 
number  of  artificial  speed  steps  and  the  spindle  speed  at  each  step  (Eq. 
13). 

CUTTING  SPEED  AND  FEED  RATE 

From  the  data  for  the  possible  number  of  feed  settings  and  the  feed 
rate  at  each  setting,  the  program  is  then  able  to  compute  Vmin  and  Vmax 
and  the  set  of  stepped  cutting  speeds  for  this  operation  (Eqs.  14,  15, 
16). 


Next,  the  constraints  on  the  feed  rate  are  established.  If  the  user 
elects  to  limit  the  range  of  feeds  suitable  for  this  cutting  operation, 
he  simply  inputs  with  the  initial  data  the  lowest  and  highest  feeds  he 

wishes  to  investigate,  as  FLOW  and  FHIGH.  But  if  the  range  is  to  include 

all  available  feeds,  FLOW  and  FHIGH  need  not  be  specified  in  the  input 
data.  In  that  case,  the  program  takes  the  first  feed  read  in  to  be  FLOW. 

The  feeds  must  be  ordered  from  low  to  high  on  the  input  data  cards. 

Constraints  on  the  cutting  speed  are  established  in  exactly  the  same 
way  as  the  constraints  on  the  feed  rate. 

TOOL  LIFE  PARAMETERS 

As  soon  as  the  initial  tool  life  data  (card  9)  are  read  into  the  com¬ 
puter  program  as  the  next  step,  all  data  supplied  or  calculated  up  to  this 
point  are  checked  for  validity.  If  a  discrepancy  is  detected,  an  error 
message  is  printed  and  the  program  is  stopped  until  the  defect  has  been 
remedied. 

When  the  program  proceeds,  NDATA  will  be  zero  on  the  first  run  and 
will  be  incremented  by  1  on  each  succeeding  run.  As  the  shop  test  data 
are  received  (card  10) ,  the  program  systematically  improves  the  tool  life 
equation  by  determining  the  number  of  data  points  available  and  applying 
the  appropriate  mathematical  steps  to  obtain  one  more  calculated  parameter 
or  to  improve  all  of  the  parameters  ALPHA,  BETA,  and  C  (see  Section  3.2, 
Steps  1  through  4) . 

When  the  value  computed  for  each  of  the  tool  life  parameters  is  equal 
to  or  less  than  the  values  specified  by  PERALF,  PERBET,  and  PERC,  the  pro¬ 
gram  accepts  the  tool  life  equation.  Subroutine  OPTM  is  then  called,  and 
the  refined  tool  life  equation  is  used  in  generating  the  final  output  of 
the  POM  computer  analysis. 

SUBROUTINE  RANK 

For  any  given  set  of  tool  life  parameters,  this  subroutine  provides 
a  ranked  set  of  optimum  feeds  and  speeds  for  the  production  objective 
specified  as  ICRIT  in  the  main  program  (card  9).  ICRIT  offers  three 
choices,  as  shown  in  the  definition  in  the  code  list  at  the  end  of  this 
section.  By  exhaustive  enumeration,  RANK  searches  all  possible  values  of 
feed,  cutting  speed,  and  depth  of  cut  to  find  those  that  best  satisfy  the 
optimization  criterion.  It  takes  into  account  all  constraints  read  into 
the  program  for  the  job,  using  only  allowable  values.  The  mathematical 
relations  for  these  calculations  are  given  at  the  end  of  Section  3.2. 
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The  output  of  RANK  each  time  it  is  called  is  two  printed  tables,  one 
showing  the  range  of  feeds  for  roughing,  the  other  the  range  for  finish¬ 
ing.  Each  of  these  tables  lists  only  the  20  best  optimizing  Values  of 
cutting  speed,  feed,  and  depth  of  cut.  From  one  or  the  other  of  the  ta¬ 
bles,  depending  on  his  purpose,  the  user  selects  the  values  of  feed  and 
cutting  speed  for  the  next  step  of  the  analysis.  He  would  logically 
choose  the  set  of  F  and  V  values  ranked  first,  but  ONLY  IF  NEITHER  HAS 
BEEN  TESTED  PREVIOUSLY.  If  one  or  both  of  those  values  occurred  in  an 
earlier  data  set,  he  will  search  down  the  list  for  a  set  of  distinctly 
new  values  and  choose  that  set. 

SUBROUTINE  OPTM 

This  subroutine  operates  in  the  same  manner  as  RANK,  but  it  is  called 
only  when  the  acceptable  values  of  tool  life  parameters  ALPHA,  BETA,  and 
C  have  been  found.  As  the  general  POM  program  is  written  (Appendix  B) , 
the  table  printouts  from  OPTM  give  the  optimum  values  of  feed,  cutting 
speed,  and  depth  of  cut  for  all  three  of  the  production  objectives  possi¬ 
ble  under  ICRIT.  Again,  there  are  two  sets  of  ranked  tables,  one  for 
roughing  conditions  and  one  for  finishing  conditions. 

Though  a  minor  amount  of  computer  time  can  be  saved  by  running  OPTM 
with  only  a  single  optimization  criterion,  the  long-range  advantage  of 
obtaining  all  of  the  information  quickly  and  simply  in  the  first  analysis 
of  a  machining  job  should  not  be  overlooked.  The  OPTM  printouts  take 
little  filing  space,  and  the  necessary  information  will  be  already  at  hand 
if  the  same  operation  should  recur^with  a  change  of  production  objective. 
This  versatility  is  particularly  desirable  when  new  materials  and  new 
cutting  tools  are  being  analyzed. 

INPUT  DATA  CODE 

The  following  tabulation  lists  and  defines  the  code  used  for  input 
data  in  the  POM  program,  with  the  corresponding  mathematical  notation. 


Input  Code 

Notation 

Definition 

NPR 

Number  or  code  to  identify  project. 

IDENT 

Type  of  lathe: 

=  1  if  speed  is  stepless.  Computer  calcu¬ 
lates  artificial  steps. 

=  0  if  speed  is  stepped.  Steps  are  read  in 
on  data  cards. 

DCMIN 

DC  . 
mm 

Minimum  depth  of  cut,  in. 

DCMAX 

DC 

max 

Maximum  depth  of  cut,  in. 

VLOW 

V. 

low 

Lowest  cutting  speed,  fpm. 

VHIGH 

Vhigh 

Highest  cutting  speed,  fpm. 

FLOW 

Flow 

Lowest  feed,  ipr. 
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Input  Code 

Notation 

Definition 

FHIGH 

Fhigh 

Highest  feed,  ipr. 

XL 

XL 

Length  of  workpiece,  in. 

DIAS 

D 

s 

Initial  diameter  of  workpiece,  in. 

DIAF 

D£ 

Final  diameter  of  workpiece,  in. 

FR 

F 

r 

Return  speed  of  carriage,  in. /min. 

Tl 

Ti 

Time  for  forward  or  backward  motion  at 
start  or  end  of  cut,  min. 

TL 

tl 

Time  for  loading  or  unloading  workpiece,  min. 

TCT 

Tct 

Tool  changing  time,  min. 

TGC 

T 

gc 

Gear  changing  time,  min. 

PRICE 

SP 

Selling  price,  $/piece. 

COSTMT 

Cmt 

Material  cost,  $/piece. 

FEEDRF 

frf 

Lowest  feed  rate  for  roughing,  ipr. 

CO 

C 

o 

Labor  and  overhead  rate,  $/min. 

CE 

C 

e 

Tool  cost,  $/edge. 

SHP 

su 

hp 

Specific  horsepower,  hp/in.  /min. 

HPM 

HP 

m 

Maximum  horsepower. 

EFF 

eff 

Efficiency  of  lathe,  decimal  form  of  percent. 

[Supply  the  following 

three  numbers  only  if  IDENT  =1.] 

RPMI 

RPMi 

Lowest  possible  speed  of  lathe,  rpm. 

RPMA 

RPM 

a 

Highest  possible  speed  of  lathe,  rpm. 

RPMCH 

RPM  . 
ch 

Desired  size  of  artificial  steps,  rpm. 

[The 

following  two  variables  apply  only  if  IDENT  =0.] 

L 

L 

Number  of  speed  steps. 

RPM(I) 

RPMi 

Values  of  L  steps,  rpm.  Must  be  ordered 
from  low  to  high. 

M 

M 

Number  of  feed  steps. 

FEED (I) 

Fi 

Values  of  M  steps,  ipr.  Must  be  ordered 
from  low  to  high. 

(CONTINUED) 
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Input  Code 

Notation 

Definition 

NDATA 

Code  for  number  of  data  points  available. 

NDATA  =  0  for  first  run  of  program. 

ICRIT 

Code  for  desired  optimization  criterion: 

1  =  minimum  unit  cost. 

2  =  maximum  production  rate. 

3  =  maximum  profit  rate. 

NVTYPE 

Code  telling  how  values  of  cutting  speed  VEL 
will  be  read  in: 

1  =  cutting  speed  in  fpm. 

2  =  cutting  speed  in  rpm. 

ALPHA 

a 

Estimated  or  previously  calculated  value 
of  tool  life  exponent. 

BETA 

3 

Estimated  or  previously  calculated  value  of 
feed  exponent . 

C 

c 

Estimated  or  previously  calculated  value  of 
tool  life  constant. 

PERALF 

ach 

Stopping  criterion  for  a.  Accept  a  <L  ac^* 

PERBET 

ech 

Stopping  criterion  for  3.  Accept  3  <.  3 

PERC 

Cch 

Stopping  criterion  for  C.  Accept  C  <.  C^. 

VEL(I) 

Cutting  speed  or  spindle  speed  tested  in  shop. 
See  NVTYPE. 

F(I) 

Feed  rate  tested  in  shop. 

T  (I) 

Tool  life  determined  in  shop  for  given 

speed  and  feed. 

3.4  PRODUCTION  EXAMPLE 

This  section  demonstrates  the  procedure  and  results  for  an  actual 
application  of  the  production  optimization  method  of  analyzing  optimum 
cutting  conditions.  Some  sample  tool  life  data  obtained  in  the  Machina- 
bility  Laboratory  at  Penn  State  were  used  for  computer  simulation  of 
production  conditions. 

PRODUCTION  INFORMATION 

The  following  initial  information  was  gathered  for  the  production 
example,  based  on  the  given  job,  the  available  tooling,  and  production 
experience  in  this  "machine  shop." 


46 


Operation  and  Tooling.  These  were  the  job  and  tool  specifications: 


Straight  turning 
Engine  lathe  (IDENT  =1) 

INCONEL-718  (Bhn-277),  24  in.  long 
C-2  carbide  insert 
8  in. 

7.5  in. 

The  tool  geometry  of  the  cemented  carbide  insert  selected  for  this  opera¬ 
tion  was  5-5-7-7-6-0-1-/16,  U.S.  standard. 

Time  Factors.  From  experience,  the  machining  time  factors  were 
established  as: 


Operation 
Machine  tool 
Workpiece 
Cutting  tool 
Diameter  at  start 
Diameter  at  finish 


Return  speed  of  carriage 
Machine  handling  time 
Work  handling  time 
Tool  changing  time 
Gear  changing  time 
for  different  rpm 


5  in. /min 
0.5  min 
5  min 
1  min 
0.5  min 


Cost  Factors.  The  labor  and  overhead  and  tool  cost  figures  were 
based  on  averages  for  this  university  machine  shop.  No  attempt  was  made 
to  approximate  industrial  figures,  which  may  vary  considerably.  The  ma¬ 
terial  cost  and  selling  price  are  assumptions  based  on  current  market 
values  and  practices. 


Labor  and  overhead  rate  $0.20/min 
Cutting  tool  cost  $0. 50/edge 

Material  cost  $300/piece 

Selling  price  $600/piece 

Machining  Constraints.  The  constraints  were  determined,  of  course, 
by  the  lathe  characteristics  and  the  nature  of  the  workpiece  and  cutting 
tool  materials. 


Maximum  horsepower  7.5 
Specific  horsepower  1.5/in. 3 
Machine  efficiency  60% 

Available  spindle  speed  range 
Available  feed  range 

Usable  cutting  speed 
Usable  feed:  finishing 
roughing 

Depth  of  cut 


/min 


20  to  1000  rpm  in  50  steps 
0.0011  to  0.0168  in  24  steps 

80.0  to  200.0  fpm 
0.003  to  0.008  ipr 
0.008  to  0.0168  ipr 

0.05  to  0.25  in. 


Tool  Life  Parameters.  For  this  demonstration  we  assumed  that  no 
data  existed  for  the  parameters  of  the  Taylor  tool  life  equation  as  it 
is  used  in  POM  analysis,  VTaF^  =  C  (Eq.  17).  The  following  figures  were 
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decided  upon  as  the  best  initial  guess,  to  be  refined  by  the  shop  tests: 

Tool  life  exponent  ct1  »  0.40 

Feed  exponent  $1  ®  0.40 

Tool  life  constant  Cj  =  15.0 

It  was  also  decided  that  a  reasonable  error  allowance  (PERALPH, 
PERBET,  PERC)  for  final  acceptance  of  the  tool  life  equation  parameters 
would  be  5%. 

Other  Specifications.  Initial  decisions  had  to  be  made  on  two  other 
points.  It  was  necessary  to  specify  the  production  objective,  and  the 
choice  for  this  case  was  minimum  cost  per  piece  (ICRIT  =>  1)  . 

Since  cutting  speed  (fpm)  can  also  be  given  as  spindle  speed  (rpm) , 
we  had  to  select  the  unit  to  be  used  fer  cutting  speed  values  from  the 
shop  tests.  We  chose  FPM  (NVTYPE  =  1). 

INITIAL  DATA  INPUT 

We  now  had  all  of  the  figures  needed  for  initial  data  input  to  the 
POM  program.  The  following  data  were  read  into  the  card  deck  in  the  order 
prescribed  by  the  format  guide  (Appendix  B) .  The  code  terms  are  defined 
at  the  end  of  Section  3.3. 


I  DENT 

= 

1 

SHP 

= 

1.5 

DCMIN 

= 

0 .05 

HPM 

= 

7.5 

DCMAX 

= 

0.25 

EFF 

= 

0.60 

VLOW 

= 

80.0 

RPM  I 

= 

20.0 

VHIGH 

= 

200.0 

RPMA 

SB 

1000. 

.0 

FLOW 

= 

0.003 

RPMCH 

= 

20.0 

FHIGH 

= 

0.0168 

M 

SB 

24 

XL 

= 

24.0 

FEED (I) 

= 

0.0011 

DIAS 

= 

8.0 

to 

O.OI681 

24  steps 

DIAF 

SB 

7.5 

NDATA 

SB 

0  (Initial  run) 

FR 

= 

5.0 

ICRIT 

= 

1 

T1 

SB 

0.5 

NVTYPE 

SB 

1 

TL 

= 

5.0 

ALPHA 

- 

0.4  ' 

TCT 

= 

1.0 

BETA 

= 

0.4 

Initial  run 

TGC 

SB 

0.5 

C 

BE 

15.0 

PRICE 

SB 

600.0 

PERALPH 

SB 

0.05 

COSTMT 

= 

300.0 

PERBET 

= 

0.05 

FEEDRF 

SB 

0.008 

PERC 

SB 

0.05 

CO 

SB 

0.2 

CE 

SB 

0.5 
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COMPUTER  PRINTOUTS 


Completing  the  optimization  of  this  example  required  five  steps  of 
computer  analysis  and  four  periods  of  shop  tests.  The  printouts  from  each 
of  the  first  four  runs  of  the  computer  program  consisted  of  three  sec¬ 
tions.  The  first  section  gave  all  of  the  job  constants,  which  were  the 
same  for  each  run,  plus  the  tool  life  conditions  pertaining  to  the  par¬ 
ticular  run.  For  the  specified  production  objective,  minimum  cost  per 
piece,  the  second  section  ranked  the  20  best  combinations  of  parameters 
for  finishing  conditions.  The  third  section  ranked  the  20  best  combi¬ 
nations  for  roughing  conditions. 

The  printout  for  the  fifth  and  final  computer  run  consisted  of  seven 
sections.  The  first  section  repeated  the  specified  conditions  and  all  of 
the  accumulated  tool  life  data,  and  gave  the  final  tool  life  parameters 
and  the  basis  for  accepting  them.  Then  followed  two  sections  for  each  of 
the  three  possible  production  objectives,  one  ranking  the  20  best  final 
optimizations  for  finishing  conditions,  the  other  the  20  best  final  opti¬ 
mizations  for  roughing  conditions. 

For  reasons  of  economy  in  producing  this  publication,  we  show  in 
Appendix  B  only  three  samples  from  the  printouts,  selected  to  indicate 
their  general  content  and  arrangement.  Reproduced  there  are  the  first 
section  and  the  finishing  table  from  Step  2,  and  the  final  optimization 
for  minimum  cost  per  piece,  finishing,  from  Step  5.  A  glance  at  the  sam¬ 
ples  at  this  point  will  give  the  reader  a  concrete  picture  of  the  form  and 
scope  of  the  computer  output,  which  will  be  useful  as  we  proceed. 

ANALYTIC  PROCEDURE 

Steps  1  through  4  in  the  following  discussion  correspond  to  the  num¬ 
bered  steps  under  PARAMETERS  OF  TOOL  LIFE  EQUATION  in  Section  3.2.  Step  5 
is  the  final  optimization. 

Step  1:  No  Shop  Data.  The  estimated  tool  life  equation  for  the 
first  run  of  the  computer  program  is 

VT°*V’4  =  15>Q 

With  that  equation,  the  first  output  of  the  program  is  the  two  tables  for 
finishing  and  roughing.  In  this  discussion  we  will  follow  the  optimiza¬ 
tion  for  only  the  finishing  condition.  Since  we  have  no  reason  to  do 
otherwise,  we  select  from  the  table  for  that  condition  the  first  ranking 
combination  to  test  on  our  shop  lathe  with  the  given  material  and  cutting 
tool : 

V  x  =85  fpm 

Fj  *  0.0078  ipr 

Dj  =  0.250  in.  (one  pass) 

As  the  average  of  several  tests  employing  that  combination  without 
change,  the  machinist  returns  to  the  program  a  tool  life  value  of 

Tj  =  26.3  min 
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This  value  can  be  determined  only  by  the  operator  in  the  production  shop, 
since  he  determines  the  time  for  tool  changes. 

Step  2:  One  Set  of  Shop  Data.  We  now  have  NDATA  =  1  (card  9),  the 
first  set  of  shop  data  V^F^Ti  to  enter  into  the  program  (card  10). 

Since  one  set  of  shop  data  permits  us  to  compute  only  one  of  the  tool 
life  parameters,  we  decide  to  hold  ALPHA  and  BETA  at  0.4  and  find  a  new 
value  C2  for  the  tool  life  equation  constant,  using  Tx  «=  26.3  min.  The 
result  of  this  computation  is 

C2  =  45.1 

and  the  tool  life  equation  becomes 
VT°  •  ‘‘F0  * **  *  45.1 

With  this  equation,  subroutine  RANK  produces  another  pair  of  tables. 
From  the  table  for  finishing  (shown  in  Appendix  B)  we  again  select  a  com¬ 
bination  to  be  tested  in  the  shop.  But  we  must  remember  that  the  values 
for  both  V  and  F  must  be  different  from  the  values  previously  tested 
(Section  3.2,  Step  2).  The  first  three  feed  values  in  the  table  are  all 
F  =  0.0078.  The  fourth  combination  shows  no  duplication  of  either  V  or 
F,  and  therefore  seems  suitable  for  testing.  So  the  combination  sent  to 
the  shop  is 

V2  =  136  fpm 

F2  =  0.0068  ipr 

D2  =  0.25  in.  (one  pass) 

The  average  figure  returned  to  the  program  from  this  group  of  tests 
on  the  lathe  is 

T2  =  7.0  min 

Step  3:  Two  Sets  of  Shop  Data.  When  NDATA  =  2,  we  can  compute  two 
of  the  tool  life  equation  parameters.  Still  holding  ALPHA  at  0.4,  we 
obtain  these  values  for  BETA  and  C: 

83  =  -0.44 
C3  =  2597.6 

and  the  tool  life  equation  becomes 
vto.4f-o.4-  =  2597  0  6 

Compared  to  the  value  computed  for  C2 ,  the  new  value  for  C3  looks  rather 
unreasonable.  But  let  us  see  what  the  program  will  do  with  it. 

Scanning  the  finishing  table  computed  with  this  equation,  we  find 
that  the  first  combination  showing  unused  values  of  both  V  and  F  is 

V3  =  14 1'.  0  fpm 

F3  =  0.0060  ipr  D3  =  0.25  in.  (one  pass) 
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With  that  set  of  cutting  conditions,  the  machinist  reports  an  average 
tool  life  value  of 

13  =  6.7  min 

Step  4:  Three  Sets  of  Shop  Data.  Now  that  NDATA  =  3,  we  can  compute 
all  three  of  the  tool  life  equation  parameters,  ALPHA,  BETA,  and  C,  and 
the  computer  gives  us 

a4  *  0.30 
$4  =  0.53 
C4  =  17.3 

The  first  thing  we  notice  about  these  figures  is  that  the  normal  operation 
of  the  program  has  automatically  brought  the  value  of  the  constant  C  back 
into  the  ballpark. 

With  the  new  tool  life  equation 
vto.30fo.53  =  17>3 

we  get  two  more  optimizing  tables.  But  this  time  our  choice  from  the 
finishing  table  does  not  necessarily  have  to  be  governed  by  V  4  Vis2,3 
and  F  ^  Fi}2,3*  Because  the  second  ranked  set  of  conditions  offers  gen¬ 
erally  better  production  values  and  only  a  slightly  higher  cost  per  piece 
than  the  set  ranked  first  in  terms  of  the  production  objective,  we  choose 
the  second  set  for  the  next  group  of  shop  tests: 

V4  =  126.0  fpm 

F4  =  0.0078  ipr 

D4  =  0.25  in.  (one  pass) 

This  combination  yields 

T4  =  7.1  min 

Step  5:  Acceptance  and  Optimization.  Calculating  that  set  of 
values,  the  computer  finds  the  tool  life  equation  parameters  to  be 

a5  =  0.30 
85  =  0.53 
C5  =  17.3 

which  means  that  the  tool  life  equation  agrees  exactly  with  the  equation 
obtained  ir.  Step  4: 

vto.3fo.53  =  17  3 

The  first  section  of  the  printout  tells  us  that  the  5%  change  limit 
specified  for  PERALF,  PERBET,  and  PERC  has  been  satisfied  and  the  program 
is  ready  to  call  0PTM,  the  final  optimizing  procedure  for  all  three  of  the 
possible  production  objectives. 
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RESULTS  AND  CONCLUSION 


Taking  the  first-rank  figures  from  each  of  the  six  tables  produced 
by  OPTM,  the  results  of  POM  analysis  of  optimum  cutting  conditions  for 
this  example  can  be  briefly  summarized  as  follows: 


Param¬ 

Min.  Unit 

Cost 

Max.  Prod 

.  Rate 

Max.  Profit  Rate 

eter 

Finishing 

Roughing 

Finishing 

Roughing 

Finishing 

Roughing 

V,  fpm 

116 

85 

126 

107 

126 

85 

F,  ipr 

0.0078 

0.0112 

0.0078 

0.0092 

0.0078 

0.0112 

D,  in. 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

T,  min 

9.25 

13.9 

7.12 

8.98 

7.12 

13.9 

That  is  a  great  deal  of  information  to  obtain  from  a  little  over  one 
minute  of  computer  time.  Note  also  that  the  20  choices  provided  by  the 
complete  tables  for  roughing  and  finishing  offer  the  industrial  user  max¬ 
imum  flexibility. 

In  this  case,  as  in  the  case  of  PIM,  we  tested  an  example  for  which 
the  optimum  cutting  conditions  had  been  determined  by  conventional  means. 
Provided  that  the  tool  life  data  returned  from  the  shop  at  each  step  are 
the  average  from  several  tests  and  the  specifications  for  choosing  the 
combination  to  be  tested  are  observed,  the  results  of  POM  analysis  are 
quite  acceptable  for  industrial  application. 
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This  symbol  denotes 


(CONTINUED) 


program 


FORMAT  GUIDE  FOR  PIM  INPUT  DATA  CARDS 


Card  Column  Number  Card  Column  Number 

Type  Variable  Number  Format _ Type  Variable  Number  Format 


7  M  1-3  13 


1 

IDENT 

1 

11 

VLOW 

11-20 

F10.0 

VHIGH 

21-30 

F10.0 

FLOW 

31-40 

F10.0 

FHIGH 

41-50 

F10.0 

VFLOW 

51-60 

F10.0 

VFHIGH 

61-70 

F10.0 

2 

SHP 

F10.0 

HPM 

11-20 

F10.0 

EFF 

21-30 

F10.0 

DC 

31-40 

F10.0 

3 

FLO 

1-10 

F10.0 

TLC 

11-20 

F10.0 

Q 

21-30 

F10.0 

D 

31-40 

F10.0 

4 

VI 

1-10 

F10.0 

VA 

11-20 

F10.0 

VCH 

21-30 

F10.0 

Include  only  if 

IDENT  = 

1. 

5 

L 

1-3 

13 

Include  only  if 

IDENT  = 

0. 

6 

L  values 

1-10 

F10.0 

of  V(I) 

to 

71-80 

F10.0 

Must 

be  ordered 

from  low 

to  high. 

If  L 

is  greater 

than  8, 

continue 

on  another  card,  same  columns. 


8  M  values  1-10  F10.0 

of  F (I)  to 

71-80  F10.0 


Must  be  ordered  from  low  to  high. 
If  M  is  greater  than  8,  continue 
on  another  card,  same  columns. 


9 

KT 

1-2 

12 

IR 

3-4 

12 

IF, NPT 

5-6 

12 

NOP ,VMAX 

11-20 

F10.0 

FOP ,FMAX 

21-30 

F10.0 

Blank 

card  on  initial  run 

(all 

values 

;  0). 

10 

VL  (I) 

1-10 

F10.0 

FD(I) 

11-20 

F10.0 

P(I) 

21-30 

F10.0 

T  (I) 

31-40 

F10.0 

TP  (I) 

41-50 

F10.0 

NPT  or 

IF  values. 

Include 

only 

when  IP 

'  ^  0.  When  IF  is  greater 

than  1 , 

continue 

on  other 

cards , 

same  columns . 


Format  15:  Integer-valued  number, 
no  more  than  5  digits  or  4  digits 
and  minus  sign,  no  decimal  point. 
Right  justified.  Read  blank  as  0. 
Examples : 


0 

1 

2 

0 

0 

read 

1200 

- 

1 

2 

0 

0 

read 

-1200 

2 

0 

0 

read 

200 

2 

read 

200 

2 

read 

2 

read 

0 

Format  F10.0;  Real-valued  number, 
no  more  than  9  digits  and  decimal 
point.  Minus  sign  may  precede. 
Computer  prints  decimal  point  in 
10th  column  if  number  has  none. 
Read  blank  as  0.  Examples: 


1 

0 

• 

1 

0 

• 

r—\ 

3 

• 

7 

3 

3 

1 

0 

0 

• 

read  10. 
read  10 . 
read  13.733 
read  100. 
read  0 
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PIM  COMPUTER  PROGRAM  LISTING 


3  «~ 

iJ  •  V 
V.  •  • 


1  z-i  *7 


cm  cm  — 

•  ♦  O 

m  ir  z 

u.  U.  o 


V  v  •  v 
*>  * 
>  T  U.  M 

»  a 

v  v  *  yj 


—  —  u.  ar 

v  S 

U-  u.  U 


—  —  CL 

I  I  o 

t  t  u. 


UJUJ^Z 

o  ^  U-  O 

<  <T  X  ►— 

tt  a  o  ►* 

^  o 
*z 
x  x  -  o 
a  a  •.  o 

u.  a  uj 

—  —  o  o 
o  o  z  z 

IK  lU  «f  - 

uj  Ui  a  w 

a.  a  vo 

l/M/)  «  uj 

cc  ►- 
n  uj  a 
Z  J  p<  j 

—  a  w  < 

*-  z  a  ~ 

K  M  UJ  K 

Da  ujm 

u  ^  u.  Z 


UJ  *  X 

a.  v  —  m 

wo  —  oc  * 


—  -  X  p- 

a  x  m  » 

VO  +  —  * 


-XXX  *> 
»  CC  00  X 

I-*  m 


—  x  *  —  * 

-  O  IX  x 
In  •  a  n 

*  -  —  a  * 

?  o  —  ^ 


•  -  •  a 

z  o 

•  •  U. 
M  f*>  * 


O'  «  <  <  <  I 

*  x  x  x  x 

-  a  a  a  a: 

♦  o  o  o  o 

•  U.  U.  U.  U.  1 

in  O  ^  o 
JA  ^  ^  N 
O'  O'  O'  O' 


•  W  X  X  M 
xw-Cp**- 
m  <  —  < 


k  x  o  a. 

V)  K  UJ  O 

a 

uj  uj  O 
X  •  o  2 

h  M 

uj  in  y 

o  a  z  ~ 
uj  uj  a  x 
►  wu  u 
u.  z  < 

VO  UJ  X 

•—  in  cc 

»-  •-  z 

<  Z  UJ 
irtO  <  > 

z  o  •- 

I/)  <-  o 
<►-»/> 

X  wo  tt  uj 

UJ  UJ  X 


o  cl  x  a 

CD  UJ  <  O 
<  X  Of  u. 
»-  < 

v  ►-  a  a  z 

-  ZD  O 

•  —  u.  to  — 

#  o  z  y- 

-  a  o  ~  — 
w  z  z  o 
^  z  —  z 

in  3  TO 

*  x  poo 

X  PM  <  < 

(O  ►*  M  I  X 

pm  a  a  D 

•  o  uj  uj  X 

-  i-  >  — 
I  UJmO^ 
*  X  cc  cd  a 

o  <r  O 
k  »  »  #  # 


5  TT  Oh< 
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[  UJ  CL  3 
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J  h  o  o  o 
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CHECK  THE  STATUS  OF  M  AND  REARRANGE  B*S  ACCOROENG  TO 
PREVIOUS  PIVOT  IF  THE  MATRIX  EQUATION  IS  TO  BE  SOLVED 


INITIAL  TESTING  CONDITIONS  RECOMMENDED: 


SP_ 

NDLE  SPEED  (N) 

FEED  (F) 

(rpm) 

TTprJ 

(1) 

160.0 

0.0068 

(2) 

220.0 

0.0068 

(3) 

260.0 

0.0068 

(4) 

160.0 

0.0102 

(3) 

160.0 

0.0092 

(6) 

220.0 

0.0102 

(7) 

220.0 

0.0092 

(8) 

260.0 

0.0102 

(9) 

260.0 

0.0092 

PROCEDURE: 

(A)  TEST  ALL  9  CONDITIONS  OF  THE  ABOVE  RECOMMENDED  COMBINATION 
OF  SPEEDS  AND  FEEDS  AT  THE  WORKSHOP  UNDER  THE  PRODUCTION 
CONDITIONS. 

(B)  COLLECT  AND  RECORD  THE  FOLLOWING  PRODUCTION  DATA  OBTAINED 
FROM  THE  SHOP  AFTER  THE  PRODUCTION  TESTS  AS  SPECIFIED: 

(1)  NUMBER  OF  PARTS  PRODUCED  (NP) 

DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(2)  NUMBER  OF  TOOL  CHANGES  (NTC)  MADE 
DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(C)  FEED  BACK  THE  ABOVE  INFORMATION  BY  PUNCHING  INPUT  CARDS 
AS  SPECIFIED  IN  THE  WRITE-UP, 

CODES  FOR  NEXT  RUN  OF  THE  COMPUTER  PROGRAM: 

IR  =  0 
IF  -  9 
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FIRST  OPTIMIZATION  ANALYSIS  BY  COMPUTER  WITH 
THE  9  INITIAL  TEST  CONDITIONS: 


THE  PERFORMANCE  INDEX  IS  DEFINED  AS: 

P I  =Q*PR+  ( 1 ..  -Q)  *  ( 1 .  /Cll)  WHERE  Q=0.0 

SPINDLE  SPEED  (N)  FEED  (F)  COST  INDEX  PRODUCTION  RATE  INDEX 


(rpm) 

( 1  PR) 

(CU) 

(PR) 

(1) 

160.0 

0.00680 

2.200 

0.05 

(2) 

220.0 

0.00680 

1.693 

0.06 

(3) 

260.0 

0.00680 

1.511 

0.07 

(4) 

160.0 

0.01020 

1.525 

0.07 

(5) 

160.0 

0.00920 

1.627 

0.06 

(6) 

220  ..0 

0.01020 

1.182 

0.09 

(7) 

220.0 

0.00920 

1.290 

0.08 

(8) 

260.0 

0.01020 

1.059 

0.11 

(9) 

260.0 

0.00920 

1.175 

0.10 

BASED 

UPON  THE  ABOVE 

PRODUCTION  TEST 

RESULTS,  THE 

FOLLOWING 

CONDITION  IS  RECOMMENDED  AS  THE  OPTIMUM  CONDITION  TO  BE  TESTED  FURTHER: 

***********  *********************** 

|  SPINDLE  SPEED  (NOP)  :  320.0  RPM  | 

*  FEED  (FOP) :  0.0094  I  PR  * 

|  CUTTING  SPEED  (VOP)  :  502.7  FPM  * 
********************************** 

FURTHER  TESTING  CONDITIONS  RECOMMENDED  BASED  UPON  THE  OPTIMUM 
CONDITION  FROM  THE  PREVIOUS  TEST  RESULTS: 

SPINDLE  SPEED  (N)  FEED  (F) 


(RPM) 

( 1  PR) 

(1) 

320.0 

0.0094 

(2). 

280.0 

0.0094 

(3) 

240.0 

0.0094 

PROCEDURE: 

(A)  TEST  ALL  3  CONDITIONS  OF  THE  ABOVE  RECOMMENDED  COMBINATION 
OF  SPEEDS  AND  FEEDS  AT  THE  WORKSHOP  UNDER  THE  PRODUCTION 
CONDITIONS 

(B)  COLLECT  AND  RECORD  THE  FOLLOWING  PRODUCTION  DATA  OBTAINED 
FROM  THE  SHOP  AFTER  THE  PRODUCTION  TESTS  AS  SPECIFIED: 

(1)  NUMBER  OF  PARTS  PRODUCED  (NP) 

DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(2)  NUMBER  OF  TOOL  CHANGES  (NTC)  MADE 
DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(C)  FEED  BACK  THE  ABOVE  INFORMATION  BY  PUNCHING  INPUT  CARDS 
AS  SPECIFIED  IN  THE  WRITE-UP. 

CODES  FOR  NEXT  RUN  OF  THE  COMPUTER  PROGRAM: 

IR  =  1 

IF  =  3  (CONTINUED) 
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FURTHER  OPTIMIZATION  ANALYSIS  BY  COMPUTER  WITH  THE  PRODUCTION 
TEST  RESULTS  UNDER  THE  PREVIOUSLY  RECOMMENDED  TEST  CONDITIONS: 

ANALYSIS  §  1 

THE  PERFORMANCE  INDEX  IS  DEFINED  AS: 


P l=Q*PR+ (l.- 

q)*(i./cu) 

WHERE  Q=0 .0 

SPINDLE  SPEED 

(N)  FEED  (F) 

COST  INDEX 

PRODUCTION 

(RPM) 

(IPR) 

(CU) 

TprT 

(1) 

320.0 

0.00940 

1.193 

0.13 

(2) 

280.0 

0.00940 

1.162 

0.11 

(3) 

240.0 

0.00940 

1.165 

0.10 

BASED  UPON  THE  ABOVE  PRODUCTION  TEST  RESULTS,  THE  FOLLOWING  CUTTING 
CONDITION  IS  RECOMMENDED  AS  THE  OPTIMUM  C®NilTI0N  TO  BE  TESTED  FURTHER: 

************************** ********* 

I  SPINDLE  SPEED  (NOP) :  300.0  RPM  * 

*  FEED  (FOP)  :  0*0102  IPR  * 

|  CUTTING  SPEED  (VOP) :  471.2  FPM  | 

*********************************** 


FURTHER  TESTING  CONDITIONS  RECOMMENDED  BASED  UPON  THE  fPTIMUM 
CONDITION  FROM  THE  PREVIOUS  TEST  RESULTS: 


SPINDLE  SPEED  (N)  FEED  (F) 


(RPM) 

(IPR) 

(1) 

300.0 

0.0102 

(2) 

260.0 

0  .0102 

(3) 

320.0 

0.0094 

(4) 

220.0 

0.0102 

PROCEDURE: 

(A)  TEST  ALL  4  CONDITIONS  OF  THE  ABOVE  RECOMMENDED  COMBINATION 
OF  SPEEDS  AND  FEEDS  AT  THE  WORKSHOP  UNDER  THE  PRODUCTION 
CONDITIONS. 

(B)  COLLECT  AND  RECORD  THE  FOLLOWING  PRODUCTION  DATA  OBTAINED 
FROM  THE  SHOP  AFTER  THE  PRODUCT!®N  TESTS  AS  SPECIFIED: 


(1)  NUMBER  OF  PARTS  PRODUCED  (NP) 

DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(2)  NUMBER  OF  TOOL  CHANGES  (NTC)  MADE 
DURING  A  SPECIFIED  TIME  PERIOD  (T) 

(C)  FEED  BACK  THE  ABOVE  INFORMATION  BY  PUNCHING  INPUT  CARDS 
AS  SPECIFIED  IN  THE  WRITE-UP, 


CODES  FOR  NEXT  RUN  OF  THE  COMPUTER  PROGRAM: 

IR  =  1 
IF  =  4 
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FURTHER  OPTIMIZATION  ANALYSIS  BY  COMPUTER  WITH  THE  PRODUCTION 
TEST  RESULTS  UNDER  THE  PREVIOUSLY  RECOMMENDED  TEST  CONDITIONS 


ANALYSIS  //  2 

THE  PERFORMANCE  INDEX  IS  DEFINED  AS: 

Pt=Q*PR+(l.-Q)*(l./CU)  WHERE  Q=0.0 


SPINDLE  SPEED  (N) 

FEED  (F) 

COST  INDEX 

PRODUCTION  RATE  INDEX 

(RPM) 

(ipr) 

(cu) 

(PR) 

(1) 

300.0 

0.  01020 

1.077 

0.13 

(2) 

260.0 

0.01020 

1.059 

0.11 

(3) 

320,0 

0.00940 

1.193 

0.13 

(4) 

220.0 

0. 01020 

1.182 

0.09 

BASED 

UPON  THE  ABOVE  PRODUCTION  TEST  RESULTS, 

THE  FOLLOWING  CUTTING 

CONDITION  IS  RECOMMENDED  AS  THE  OPTIMUM  CONDITION  TO  BE  TESTED  FURTHER: 

************************************ 

|  SPINDLE  SPEED  (NOP) :  300.0  RPM  | 

*  FEED  (FOP) :  0.0102  IPR  * 

*  CUTTING  SPEED  (VOP) :  471.2  FPM  * 

********************** ************** 


******************************************************* 
|  THE  OPTIMUM  POINT  ABOVE  HAS  SATISFIED  THE  STOPPING  | 

*  CRITERION.  NO  FURTHER  TESTING  IS  NEEDED.  THE  * 

*  ABOVE  MACHINING  PARAMETERS  CAN  BE  CONSIDERED  AS  THE  * 

*  OPTIMUM  CONDITION  FOR  THE  GIVEN  MACHINING  OPERATION,* 
******************************************************* 
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FLOW  DIAGRAM,  POM  COMPUTER  PROGRAM 


< 

►-  TJ 
<  i- 
O  (Q 
Z  O 
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ICRIT 
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Also  test  for  BETA.C 


FORMAT  GUIDE  FOR  POM  INPUT  DATA  CARDS 


Card 

Column 

Number 

Card 

Column 

Number 

Type 

Variable 

Number 

Format 

Type 

Variable 

Number 

Fo rmat 

1 

NPR 

1-5 

15 

6 

L  values , 

1-5 

F5.0 

IDENT 

6-10 

15 

RPM(I) 

to 

DCMIN 

11-20 

FlOoO 

76-80 

F5.0 

DCMAX 

VL0W 

21-30 

31-40 

F10.0 

F10.0 

Must 

be  ordered 

from  low 

to  high. 

VHIGH 

FLOW 

41-50 

51-60 

F10.0 

F10.0 

7 

M 

1-3 

13 

FHIGH 

61-70 

F10„0 

8 

M  values , 

1-5 

F5o0 

2 

XL 

DIAS 

1-10 

11-20 

F10.0 

F10.0 

FEED (I) 

to 

76-80 

F5.0 

DIAF 

21-30 

F10.0 

Must 

be  ordered 

from  low 

to  high. 

FR 

31-40 

F10 ,0 

If  M 

is  greater 

than  16 , 

continue 

T1 

41-50 

F10.0 

on  another  card 

,  same  columns. 

TL 

51-60 

F10.0 

TCT 

61-70 

F10  .0 

9 

NDATA 

1-10 

110 

TGC 

71-80 

F10.0 

ICRIT 

NVTYPE 

11-20 

21-30 

110 

110 

3 

PRICE 

1-10 

F10.0 

ALPHA 

31-40 

FlOoO 

CO  STMT 

11-20 

F10.0 

BETA 

41-50 

FlOoO 

FEEDRF 

21-30 

F10.0 

C 

51-60 

F10.0 

CO 

31-40 

F10.0 

PERALF 

61-65 

F5.0 

CE 

41-50 

F10.0 

PERBET 

66-70 

F5 . 0 

SHP 

HPM 

51-60 

61-70 

F10.0 

F10.0 

PERC 

71-75 

F5.0 

EFF 

71-80 

F10..0 

10 

VEL(I) 

F(I) 

1-10 

11-20 

FlOoO 

FIO.,0 

4 

RPMI 

1-10 

F10 . 0 

1 

T  (I) 

21-30 

F10.0 

RPMA 

RPMCH 

11-20 

21-30 

F10,0 

FlOoO 

Include  only  when  NDATA  f 
One  card  for  each  set  of 

0, 

data. 

Include  only  if 

IDENT  = 

1. 

5 

L 

1-3 

13 

Include  only  if 

IDENT  = 

0. 

Format  15:  Integer-valued  number, 
no  more  than  5  digits  or  4  digits 
and  minus  sign,  no  decimal  point. 
Right  justified.  Read  blank  as  0. 
Examples : 


0 

1 

2 

0 

0 

read 

1200 

- 

1 

2 

0 

0 

read 

-1200 

2 

0 

0 

read 

200 

2 

read 

200 

2 

read 

2 

read 

0 

Format  FlOoO:  Real-valued  number, 
no  more  than  9  digits  and  decimal 
point.  Minus  sign  may  precede. 
Computer  prints  decimal  point  in 
10th  column  if  number  has  none. 
Read  blank  as  0.  Examples: 


1 

0 

1 

o 

1 

3 

• 

7 

3 

3 

1 

0 

0 

c 

read  10. 
read  10 „ 
read  13.733 
read  100, 
read  0 
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TOOL  LIFE  EXPONENT:  ALPHA  •  0.4000 
FEED  EXPONENT:  BETA  -  0.4000 
THE  TOOL  LIFE  CONSTANT,  C,  IS  ESTIMATED  TO  BE: 

C  -  45.083 
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Both  methods  are  rapid  and  economical,  and  the  programs  can  be  processed  by  Both  methods  are  rapid  and  economical,  and  the  programs  can  be  processed  by 
either  in-plant  or  remote  computer  facilities.  The  user  is  given  all  infor-  either  in-plant  or  remote  computer  facilities.  The  user  is  given  all  infor¬ 
mation  needed  to  install  the  programs  and  adapt  them  to  his  purposes.  mation  needed  to  install  the  programs  and  adapt  them  to  his  purposes. 


